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DIGITAL SIGNAL PROCESSING 
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Teaching Schedule 
Hours/Week 


Examination Scheme 


Theory 


Tutorial 


Practical 


Internal Assessment 


Final 


Total 
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3 


- 
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Theory 


Practical* 


Theory 


Practical** 


20 


25 


80 


- 



* Continuous 

** Duration: 3 hours 
Course objectives: To provide 

1. Discrete signals 

1.1 Discrete signals - unit impulse, unit step, exponential sequences 

1.2 Linearity, shift invariance, causality 

1.3 Convolution summation and discrete systems, response to discrete inputs 

1 .4 Stability sum and convergence of power series 

1.5 Sampling continuous signals spectral properties of sampled signals 

2. The discrete Fourier transforms 

2. 1 The discrete Fourier transform (DFT) derivation 

2.2 Properties of the DFT, DFT of non-periodic data 

2.3 Introduction of the fast fourier transform (FFT) 

2.4 Power spectral density using DFT/FFT algorithms 

3. Z transform 

3.1 Definition of Z transform one sided and two sided transforms 

3.2 Region of convergence relationship to causality 

3.3 Inverse Z transform - by long division, by partial fraction expansion. 

3.4 Z transform properties - delay advance, convolution, Parseval's theorem 

3.5 Z transforn transfer function H (Z) -transient and steady state sinusoidal 
response pole zero relationships, stability 

3.6 General form of the linear, shift invariant constant coefficient difference 
equation 

3.7 Z transform of difference equation. 

4. Frequency response 

4. 1 Steady state sinusoidal frequency response derived directly from the 
difference equation 

4.2 Pole zero diagrams and frequency response 

4.3 Design of a notch filter from the pole zero diagram. 



Discrete filters 6 

5.1 Discrete filters structures, second order sections ladder filters frequency response 

5.2 Digital filters finite precision implementations of discrete filters 

5.3 Scaling and noise in digital filters, finite quantized signals quantization error 
linear models. 



6. HR Filter Design 7 

6.1 Classical filter design using polynomial approximations - Butterworth Chebishev 

6.2 HR filter design by transformation matched Z transform impulse, invariant transform and 
bilinear transformation 

6.3 Application of the bilinear transformation to HR low pass discrete filter design 

6.4 Spectral transformations, high pass, band pass and notch filters. 

7. FIR Filter Design 3 

7.1 FIR filter design by fourier approximation the complex fourier series 

7.2 Gibbs phenomena in FIR filter design approximations, applications of window 
Functions to frequency response smoothing rectangular hanning Hamming and Kaiser 
windows. 

7.3 FIR filter design by the frequency sampling method 

7.4 FIR filter design using the Remez exchange algorithm 

8. Digital filter Implementation 

8.1 Implementations using special purpose DSP processors, the Texas Instruments TMS320. 

8.2 Bit serial arithmetic distributed arithmetic implementations, pipelined implementations 



Laboratory: 

1. Introduction to digital signals sampling properties, aliasing, simple digital notch filter behaviour 

2. Response of a recursive (HR) digital filter comparison to ideal unit sample and frequency 
response coefficient quantization effects. 

3. Scaling dynamic range and noise behaviour of a recursive digital filter, observation of nonlinear 
finite precision effects. 



Digital Signal processing: - 

Application: - 

1) Speech recognition. 

2) Telecommunication. 



Digital Signal processing over analog signal processing: 

1) Accuracy. 

2) Offline processing. 

3) Software control. 

4) Cheaper than analog counterpart. 

Basic element of D.S.P system: - 



Analog 
input signal 



A/D 

Converter 



Digital 

signal 

processing 
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Converter 



Analog 
output signal 



Digital 
input signal 



Digital 
output signal 



Fig: Block diagram of DSP 

Most of the signal encountering science and engineering are analog in nature i.e the signals are 

function of continuous variable substance in usually take on value in a continuous range. 

To perform the signal processing digitally, there is need for interface between the analog signal and 

digital processor. This interface is called analog to digital converter. The o/p of A/D converter is 

digital signal i.e appropriate as an i/p to the digital processor. 

Digital signal processor may be a large programmable digital computer or small microprocessor 

program to perform the desired operation on i/p signal. 

It may a also be a hardwired digital processor configure to perform a specified set of operation on 

the i/p signal. 

Programming machine provide the flexibility to change the signal processing operation through a 

change in software whereas hardwired m/c are difficult to reconfigure. 

In application where the distance o/p from digital signal processor is to be given to the user in analog 

form, we must provide another interface on the digital domain into analog domain. Such an interface 

is called D/A converter. 



Advantage of Digital over analog signal processing: 



1. 



2. 
3. 



A Digital programmable system allows flexibility in reconfiguring the digital signal processing 
operation simply by changing the program. Reconfiguration of analog system usually implies 
redesign of hardware followed by testing and verification to see that if operates properly. 
Digital system provide much better control of accuracy requirements. 

Digital system are easily stored on magnetic media without loss of signal beyond that 

introduce in A/D conversion. As a consequence, the signals become transportable and can be 
processed offline in a remote laboratory. 



4. Digital signal processing method also allows for the implementation of more sophisticated. 

5. In some cases a digital implementation of signal processing system is cheaper than its analog counter 
part. 

* Signal: 

It is defined as any physical quantity which is a function of one or more independent variable and 
contains some information. 

In electrical sense, the signal can be voltage or current. The voltage or current is a function of time 
as an independent variable. 

The independent variable in the mathematical representation of a signal may be either continues or 
discrete. 

Continues time signals are defined analog continues times. Contineous time signals are often 
referred to as analog signals. 

Discrete time signals are defined as certain time instant. 
Digital signals are those for which both time and amplitude discrete. 

Discrete- time signals or sequence: - 

D.T.S are re mathematically as a sequence of numbers. A sequence of numbers x in which 

the n th no in the sequence is denoted by x(n) and written as: 
x = { x(n) } - infinity < n < infinity 

Figure: 



unit sample or unit impulse sequence: 

It is denoted by 5(n) and defined as: 

5{n) = 



1 



n = \ 



* Unit step sequence: 

It is denoted by U(n) and is defined as: 

U(n) = l ,n>0 
= 0, n<0 



12 3 



U(n) = 5(n) + 5(n-l) + 5(n-2)+ 



= £<?(n-*) 



U(n)= f^Sik) 



*=— « 
i.e the value of unit step sequence at time n is equal to the accumulated sum of value at index n and all 
prvious value of impulse sequence. 

Conversely the impulse sequence can be expressed as the first backward difference of unit step 
sequences. 

i.e. 8(n) = u(n) - u(n-l) 

* Unit ramp sequence:- 

It is denoted by u r (n) and defined as 
u r (n) = n n > 
= n<0 



.j^ 



* Exponential Sequence:- 

The exponential signal is a sequence of the form x(n) = a n for all n. 
If the parameter 'a' is real, then x(n) is real signal. Fig illustrate x(n) for various values of parameter 



a>l 



0<a<l 



a<] 



Fig: Graphical representation of Exponential signals. 



0<a< 1: 

Eg a = Vi 

i.e (l/2) n = 0, Vi , %, 1/8 (exponential decreasing) 

(see figure (i)) 



a> 1, 

E.g,a = 2 

i.e (2) n = 1,2, 4, 8 (Exponential increasing) 

(See figure (ii)) 

-l<a<0 
Eg.: -1/2 

i.e (-l/2) n = 1,-1/2, Va, -1/8 
(Seefig(iii)) 

a<-l 

E.g a = -2 
(-2) n =l,-2,4,-8 
(See figure(iv) ) 

Exponential sequence:- 

When the parameter 'a' is complex values , it can be expressed as: 
a = re J 
Where r and are new parameters. Hence, we can express x(n) as: 
X(n) = r n e je 

= r n (cos On + jsin 0n) 
Since, x(n) is now complex values, it can represented graphically by plotting the real part, 
x e (n) = r n cos On as a function of n and separately plotting the imaginary part. 
xj(n) = rsinBn as a function of n. 
Fig. illustrates the graphs fo 



X e (n) and x ; (n) . 



I | • . ' | 

• ■ 



We observe that the signals x e (n) and xi(n) are damaged (decaying exponentially, i.e r <r ) cosine 
function and damped sine function. 
If r =1 , the damping disappears and x e (n) , x;(n) and x(n) have fixed amplitude which is unity. 

Alternatively, the signal x(n) can be represented graphically by the amplitude function. 
Ix(n)l = A(n) = r n 
And phase function 

Zx(n) = 0(n) = 9 n 
Representation of discrete-time signal: - 

1) Functional representation :- 

E.g x(n) = 1 , for n = 1,3 
= 4, for n = 2 
= elsewhere, 

2) Tabular representation: 



n 


2 


-1 





1 


2 


3 


4 


x(n) 











1 


4 


1 






3) Sequence representation: 

X(n) = { 0,0, 1,4, 1,0, } Infinite duration . 

X(n) = { 0, 1, 4, 1) finite duration (4- point sequence) 

4) Graphical representation :- 

Figure: 
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Linearity: A system is called liner of superposition principal applies to that system. This means that 
the liner system may be defined as one whose response to the sum of weighted inputs is same as the sum 
of weighted response. 

Let us consider a system. If xi(n) is the input and yi(n) is the output. Similarly y2(n) is the response to 
X2(n) . Then for liner system. 

a^in) + a l x 2 (n) — > a^^n) + a 2 y 2 (n) (1) 

For any nonlinear system the principle of superposition doesnot hold true and equation (i) is not 
satisfied. 



Numerical: 

For the following system, determine whether the system is liner or not. 
(1) y(n) = 2x(n)+3 
Solution: 

yi(n) = 2xi(n)+3 
y 2 (n) = 2x 2 (n)+3 



Applying superposition principle suppose, 

x(n) = aixi(n)+a 2 x 2 (n) 
Then, 
y 3 (n) = 2[aiXi(n)+a 2 x 2 (n)]+3 
= 2aiXi(n)+ 2a 2 x 2 (n)+3 
y 4 (n) = aiyi(n)+ a 2 y 2 (n) 

= ai [2 Xl (n)+3]+ a 2 [2x 2 (n)+3] 
= = 2xi(n)ai +2a 2 x 2 (n)+3ai+3a 2 
Since, 

y 3 (n) & y 4 (n) . The system is nonlinear. 

(ii) y(n) = x 2 (n) 
Solution: 

y 1 (n) = xf(n) 
y 2 {n) = x 2 2 {n) 
Applying input such that, 
x(n) = aixi(n)+a 2 x 2 (n) 
The response will be, 
y 3 (n) = [aixi(n)+a 2 x 2 (n)] 2 
and again, 
y 4 (n) = aiyi(n)+a 2 y 2 (n) 

= a l xf(n) + a 2 yl(n) 
.'. y 3 (n) yt y 4 (ft) The system is nonlinear. 

Shift invariance:- 

A system is shift invariant, if the input output relationship doesnot vary with shift. In other words for a 
shift invariant system shift in the input signal results in corresponding shift in output. Mathematically, 

x(n) >y(«) 

Which means that y(n) is the response for x(n). If x(n) is shifted by n , then output y(n) will also be 
shifted by same shift no i.e 

x(n -n ) > y(n -n Q ) 

Where n is an integer. 

If the system doesnot satisfy above expression, then the system is called shift variant system. 

The system shifting both linearly and time invariant properties are popularly known as liner time 
invariant system or simply LTI systems. 

Numerical: 
Check whether the system are shift invariant or not. 
(i) y{n) = x\n) 
Solution: 

Let us shift in input by n , then the output will be, 



y l (n) = x (n-n Q ) 
.:y l (n) = y(n-n Q ) 
Hence the system is shift invariant system. 

(ii) y(n) = x(mn) 

Let us shift in input by no then, response will be, 

yi(n) = x(mn - n ) 
and the shift in output by n will yields. 
y(n — n ) = x[mn - mn ] 
:. The system is shift variant. 

Causality: - 

A system is causal of the response does not begin before the input function is applied. This means that 
of input is applied at n= no , then for causal system output will depend on values of input x(n) for n<n Q . 
Mathematically, 

y(n ) = T[x(n), n<n ] (i) 

The response of causal system to an input doesnot depend on future values of that input but depends 
only on present or past values of input. 

On the other hands, of the response of the system to an input depends on future values of the input, 
the system is noncausal. A non causal system doesnot satisfy equation (i). 

Causal system are physically realizable whereas non causal system cannot be implemented practically. 
There is no system possible practically which can produce its output before input is applied. 
The equation, 
y(n) = x(n - 1) describes the causal system and 
y(n) = x(n) - x(n + 1) describe the non causal system. 



Memory less system :- 

A system is referred to as memory less if the o/p y(n) at every value of n depends only on the i/p x(n) 
at the same value of n. 
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* Response of LTI system to Arbitrary input convolution solution:- 

Consider any arbitrary discrete time signal x[n] as shown in fig: 



->" 



A discrete time sequence signal be represented by a sequence of individual impulses as shown in figure: 



x(0)5(n) 



x(l)5(n-l) 



1 2 3 



We can write, 

X(n) = +x(- 1 )6(n+ 1 ) +x(0) 5(n) +x( 1 ) 8(n- 1 ) +x(2) 5(n-2) + 

= ^xik^in-k) (i) (convolution sum) 



Suppose h(n) is the o/p of LTI system when 8(n) is i/p. Therfore, the o/p for i/p x(-l) 8(n+l) is x(-l) 
h(n+l). Then the o/p y(n) for the i/p x(n) given in equation (i) will be, 



y(n) = '^ i x(k)h(n - k) 



.(») 



k=- 



Symbolically, 

y(n) = x(n)* h(n) (iii) 

y(n) = h(n)*x(n) (iv) 



Numericals: 

* The impulse response of invalid time response is: h(n) = {1,2, 1,-1} x(n) = {1,2,3, } 
Solution: 
The response of the LTI system is , 

y(n) = ^x(k)h(n-k) 

jfc=-oo 

For n = , 



y(0)= Y^xikM-k) 



k=- 



£x(kM-k) 



k=o 



= x(0)h(0) + x(l)h(-l) 

=2+2=4 



For n = - 1 

y(-l)= f^xikM-l-k) 



k=- 



->k 



2 2 



12 3 



U 



>k 



->k 



>k 



n(-l-k) 



-2 -1 



For n = -2 

Y(-2)= £*(*)ft(-2-*) 



k=- 



x(k)n(-l-k) 



Again testing for +ve side: 

For n= 1 
Y(l)= ^x(k)h(l-k) =1+4+3 = 8 



o i 2 



-l 



For n = 2 

Y(n)= ^x(k)h(2-k) 



k=~~ 



-1+2+6+1 = 8 

3 

2 

1 



->k 



1 2 3 



■>k 



-l 



For n = 6 = 
Y(n)={ 0,1,4,8,8,3,-2,-1,0 ) 

Figure: 



x(k)(l-k) 
>k 



1 2 3 
x(k)h(l-k) 



->k 



* Determine the o/p y(n) of relaxed linear time invariant system with impulse response: 
h(n) = a" u(n) , \u\ < 1 
When the i/p is unit step sequence ie; 
x(n) = u(n) 

Solution: - 

y(n) = ^x(k)h(n-k) 
k=-°° 
For n = , 

y(0)= £*(*)/»(-*) 

k=-°= 



1111 



a2 



12 3 



■>k 



12 3 



->k 



a 2 



a 3 



h(-k) 



-3 -2 -1 



->k 



For n =1: 



y (l)= Yj^ikMl-k) =l+a 



a 2 



a2 



hd-k) 



-10 12 



"^k 



x(k)h(l-k) 



->k 



For n = 2: 



y(2) = £ x(n)/?(2 - n) = 1+a+a 2 



£=-<* 



a2 



1 2 



x(k)h(2-k) 



^k 



s= ad-r") 

1-r 
For n > 

y(n) = 1+a 2 + 

1-a" 



.+a" 



1-a 



For n < 
y(n) = 

y(oo) = lim n — > oo y(n) = lim n — > °° 



1 — a 



n+l 



1-a 



1 



1-a 




1+a+a 2 



1 2 



Interconnection of LTI system: - 



^ hl(n) ^ h2(n) >y(n) 



-3> hl(n)*h2(n) J>y(n) 



When two LTI systems with impulse response hi(n) and h 2 (n) are in cascade form the overall impulse 
response for the cascaded 2 impulse system will be, 
h(n) = hi(n) * h 2 (n) (i) 



* The parallel combinations of LTI systems and equivalent system is shown below: 



x(n)- 



> hl(n) 



> h2(n) 



+J— >yw 



> hl(n)+h2(n) >y(n) 



Determine the impulse response for the cascade of two LTI systems having impulse responses. 
hj(n) = (l/2) n u(n) 
hi{n) = (1/uf u(n) 

Solution: - 
The overall impulse response is 
h(u) = hi(n)* h 2 (n) 

= 2_ l h l (k)h 2 (n-k) 



1/2 



hl(k) 



1/4 



1/4 



h2(k) 



1/6 



->k 



->k 



1/4 



1/6 



h2(-k) 
1 



->k 



1/4 



1/6 



h(n-k) 
4 1 

n>0 



->k 



\n— & 



= £(l/2)*(l/4)" 

= (l/4) n £(l/2)*(l/4)~ 



*=() 



(l/4) n £(2) A 



= (l/4) n 






(l/4) n (2 n+1 -l) 



(l/2) n [ 2 - (1/2) 11 ] ,n>0 



Note:- If we have 'L' LTI system is cascade with impulse responses h^n) and h 2 (n) h L (n) , the 

impulse response of equivalent LTI system is 
h(n) = tn(n) * h 2 (n)*h 3 (n) *h L (n) 

Discrete system response to discrete input: - 

For a discrete-time system, consider as input sequence x(n) = e Jwn , for -go <n< qo the output of LTI 
system with impulse response h(n) is, 

y(n) = ^h(k)x(n-k) 

= ^h(.k)e M "^ k) 

k=-°° 

= e jwn Y^hiky-^ 



If we define, 
Hle^ ^/z(fcK M 



k=- 



jwi jwn 



Then, y(n) = Hie 1 "! e 

Hle Jw l describes the change in complex amplitude of complex exponential as a function of frequency w. 

H(e JW ) is called frequency response of the system. 

In general H(e JW ) is complex and can be expressed in terms of its real and imaginary parts as: 

H(e jw ) = H R (e JW )+jH I (e JW ) 
Or , In terms of magnitude and phase as, 



H(e jw ) = IH(e jw )le j ^ H(ejw) 
Let, 
x(n) = Acos(w n +0) 

= A/2 e jwon .e j0 +A/2 e - jwon . e" j0 
The response to xi(n) = A/2 e J0 e Jwon is, 

yi(n) = H(e jw )A/2e j0 e jwon 
The response to x 2 (n) = A/2e" J0 e" Jwon is , 

y 2 (n) = H(e" jw ) A/2 e" j0 . e- jwon 
Thus, total response is , 

y(n) = A/2 [H(e jwo ) e - j0 . e jwon +H(e- jwo ) e - j0 . e - jwon ] 
= A IH(e" jwo )l cos(w o n+0+9). 

Where, 9 = Z H(e jwo ) 
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Stability: We defined arbitrary relaxed system as BIBO stable. If an only if o/p sequence y(n) is 
bounded for every bounded i/p x(n). 

If x(n) is bounded their exist a constant M x such that, 
|jc(n)| <M X <oo 
Similarly if o/p is bounded their exists a constant M y such that 
\y(ri) < M < oo for all n, 
Now given such a bounded input sequence x(n) to LTI system with convolution formula. 



y(n) = 



y]h(k)x(n -k) 
< Y^\h(k)x(n-k) 

< ^\h(k)\\x(n-k)\ 



k=- 



<i^W 



k=- 



From this expression we observe that the system is bounded if the impulse response of the system 
satisfied condition, 



S„=f>«|< 



k=- 



That is linear time invariant system is stable if its impulse response is absolutely summable. 

# Determine the range of value of the parameter a for LTI system with impulse response. 

h(n) = a n u(n) is stable. 
Solution, 



S„ = £|/i(*)|<~ 

= Y\h{k)\ = Y\a k 



k=0 



k=0 



= l + \a\ + \a 2 \ + \a 3 \ + . 



a — \a\ 



Provided that \a \< 1 



Therefore the sytem is stable if a < 1 otherwise it is unstable. 

# Determien range of value of a,b for which LTI system with impulse response. 
h(n) = a" n>0 

= b n n<0 is stable. 

S n = f J \h(k)\=f J \b k \ + f j \a k \ 

k=-°* t=-~ k=0 

The first sum coverage for \a < 1 . The second sum can be manipulated as, 



k=- 



1 1 



k=\ \b\ 
1 



H \b\ b 



r 



, 1 1 



= 1 + j3 + j3 2 +j3 3 + 

= J3(1/1-J3) Provided that , 1^1 < 1 or \b\>\ 



Hence the system is stable is both \a\ < 1 or \b\>l 



# Determine whether the given system is BIBO stable or not. 

y (n) = 1/3 [ x(n) +x(n-l) +x (n-2)] 

Solution: - 
Assume that, 
lx(n)l < M n < co for all n. 

Then I y(n)l = 1/3 [ x(n) +x(n-l) +x (n-2)] 

< l/3[lx(n)l+lx(n-l)l+lx(n-2)l] 

< 1/3 [ M x + M x + M x ] 

< M x 
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BIBO = 

Bounded input 
bounded output. 



Since, M x is finite value , ly(n)l also finite. Hence the system is BIBO stable. 

# # Determine whether the given system is BIBO stable or not. 
Y(n) = r n x(n) where r>l. 

Solution: 
Assume 

\y(n)\ = \r"x(n)\ 

= r \x(n)\ 

With r >1, the multiplying factor r n diverges for increasing n and he=nce o/p not bounded. Hence the 
system is BIBO unstable. 

Nyquist Sampling theorem: - 

A signal whose spectrum is band limited to B Hz (G(w) = for Iwl > 2nB ) can be reconstructed 
exactly form its samples taken uniformly at the rate R > 2BHz (sample/sec) . In other words, minimum 
sampling frequency f s = 2Bhz. 

Consider a signal g(t) whose spectrum is band limited to Bhz. Sampling g(t) at the rate of f s hz can be 
accomplished by multiplying g(t) by impulse 5t s (t) consisting of unit impulses repeating periodically 
every T s second. Where, T s = l/f s 

Figure; 



Trigonometric fourier series of impulse train, 

S T (t) = — |l + cosw t + 2cos2w t + 2cos3w t + | 

T, Where, W s = 2ti/T s 

g(t)=g(t)S Ts (t) 

= — [g(t) + 2g(t) cos w s t + 2g(t) cos 2wJ + 2g(t) cos 3wJ + ] 

Ts 

Using modulation property, 

2g(t) cosw s t <-> (F.T) G(w-w s )+G(w+w s ) 

G(w) = — [G(w) + G(w- w s ) + G(w + w s ) + ] 

s 

= — ^G(w-nw s ) 

If we want to reconstruct g(t) from g(t) bar we should be able to recovered G(w) form G (w) . This is 

possible if there is no overlap between successive cycle of G (w). Figure (e) shows that this requires f s 
greater then 2B. 



From figure we see that g(t) can be recovered form sample g(t) by passing sampled signal through 
ideal low pass filter with bandwidth B hz. 
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Sampling of analog signals: 

There are many ways to sample analog signal. We limit our discussion to periodic or uniform sampling 
which is he types of sampling used most often in topic in practice . This is described by the relation. 
X(n) = x a (nT), - oo < n < go 

Where x(n) is discrete-time signal obtained by talking samples of analog signal x n (t) every T second. 
T = sampling period or sample interval. 

F s = 1/T = damping frequency or sampling rate (sample/sec or Hz) 
t = nT = n/F s 

In general the smapling of continuous time sinusoidal signal 
x a (t) = Acos(2nf t + 0) 

With sampling rate f s = — results discrete time signal. 

x(n) = Acos(2nft + 0) 

f 
Where, / = — - = Relative frequency of sinusoid. 

J s 

Q. Consider the analog signal x a (t) = 3cosl00nt 

a) Determine the minimum sapling rate require to avoid aliasing. 

b) Suppose that the signal is sampled at the rate F s = 200 hz What is the discrete time signal obtained 
after sampling. 

c) Suppose that the signal is sample at the rate F s = 75 hz what is the discrete time signal obtained after 
sampling. 

d) What is the frequency 0<f< FJ2 of a sinusoids that yields samples identical to those obtained in part 
c. 

Solution: - 
X k (t) = 3cosl00 7rt 
2 Ttfo = 100 % 
F = 50 hz 

b) Mininum sampling rate = 100 hz 

c) F s =200hz 

d) F s = 75hz 

X(n) = 3cos(100 n n/f s ) = 3cos(100 n n/75) = 3cos(4 jm3) 

= 3cos2 jc(2/3)n 
=3cos2 7i(l -l/3)n 
= 3cos2 7t(l/3)n 

e) F s = 75 hz , f= 1/3 
f=F /F s ,F = fF s = 25hz 



y(t) = 3cos50 Kt 

Q. Consider the analog signal at x a (t) = 3cos50 ffi + 10sin300 M-coslOO M. What is the Nequest rate 

for this signal. 

Solution: - 

Frequency present in the signal, 

Fi = 25hz, F 2 = 150 hz , F 3 = 50 hz 

F ma x= 150 hz 

Nyquest rate = 2F max = 300 hz. 



Q. Consider the analog signal at x a (t) = 3cos 2000M+5sin6000 M+10cos2000 M. 

a) What is the Nyquest rate for the signal. 

b) Assume now that we sample the signal using sampling rate F s = 5000 samples per second. What is 
the discrete time signal obtained after sapling. 

Solution: 

a) Fi = lkhz F2 = 3khz F3 = 6khz 
Nyquest rate =12 khz 

b) F s = 5000 hz = 5khz 

x(n) = 3cos(2000 n n/F s ) + 5sin(6000 n n/F s )+10cos( 12000 n n/F s ) 
= 3cos(2 7in/5) +5sin(6 Ttn/5) + 10cos(12 7tn/5) 
= 3cos(2 7in/5)+ 5sin2 7i(3/5)n+10cos2 7r(6/5)n 
=3cos(2 7m/5)+ 5sin2 7i(l-2/5)n+10cos2 7t(l+l/5)n 
=3cos(2 ko/5)- 5sin2 7i(2/5)n+10cos2 7t(l/5)n 
=13cos(27in/5)-5sin27i(2/5)n 
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Chapter:- 2 

Discrete Fourier transform :- 

Frequency domain sampling: Discrete Fourier transform. (DFT) 

Let us consider aperiodic discrete-time signal x(n) with Fourier transform 

X(w)= f>(n)<T^ (i) 

Suppose that we sample X(w) periodically in frequency at a spacing of 5(w) radian between the 
successive samples. Since X(w) is periodic with period 2n only samples in the fundamental frequency 
range are necessary. We take N equidistant sample in the interval <= w <=2n . With sample spacing 

5w = 2ti/N 




0*1 K- k w 
8(0 



Fig: Freq. domain sampling. 



If we evaluate, equation (i) we get, 
W = 2ti/N. K, k = 0toN-l 

x(—)= fx(n)e- i2 * nlN 

Summing in equation (2) can be subdivided into infinite number of summations where each sum 
contains N terms Thus, 



X 



f 2nk^ 



V N j 



2N-\ 



3N-1 



, V , \ -i2nkn/N , V / n -jlnknlN , V , \ -jlAnlN , 



n=0 



n=2N 



IN+N-X 



X ^xMe-^ 11 



l=-c° n— In 



If we change the index in the inner summation form n to n-ln and integrating the order of summation we 
obtained, 



X 



2M_ 

N 



N-\ 

= 1 



N-\ 



2]x(n-nl) 



-jlaknl N 



•0) 



The signal 

X p (n) = ^x(n-nl) 



•(4) 



/=- 



Obtained by period representation of x(n) every N samples is clearly period with fundamental period N. 

Since x p (n) is period extension of x(n) given by equation (4) it is clear that x(n) can be recovered 
from x p (n) if there is no alising in the time domain that is x(n) is limited to less than the period N of 
x P (n) 

,x(n) 



1 P 



x D (n) 



N>L 



r_ 



^_ 



L N 



r_ 



->n 



AY 



x D (n) 



N<L 



1 



1 



1 



->» 



L-l 



In summary a final duration sequence x(n) of a length L as fourier transform X(w) = ^x(n)e JH 



n=0 



for 0<=w<=2ti (5) 

When we sample X(w) at equal space frequency w k = 27iK/N , k = 0,1, 2 N-l 

Where, N > L, 

The resultant sample are, 



X(k) = X 



v 



N 



x(n)e 



(6) 



j «=o 
Where,k = 0, 1, N-l 



The relation in equation (6) is a formula for transforming sequence x(n) of length L <= N into a 
sequence of frequency samples {X(k)} of length L. Since the frequency samples are obtained by 
evaluating the fourier transform X(w) at a set of N equally spaced discrete frequencies. The relation in 
equation (6) is called discrete fourier transform (DFT) of x(n). 



x p (n) can be written as, 



N-l 



x p( n )= Yj a k ei 



2 An IN 



n = 0, 1, N-l ....(7) 



k=0 



With Fourier coefficient, 



S L k =l/N^x p (n)e- j2 '* n ' N k = 0,l, N-l 



•(8) 



k=0 



Form equation (3) and (8) 



a k = 1/N 
Therefore 



f 27^ 



k = 0, 1, N-l 



•(9) 



lV-1 



x p (n)=l/N Y,X 



Ink 
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e J2AnlN n = 0> x N _j (1Q) 

k=Q V N , 

This relation allows us to recover the sequence x(n) from frequency sample. 



x(n)= 1/N Y j X{k)e j2AnlN n = 0, 1 N-l (11) 



k=o 



This is called inverse DFT (ID FT) ( 



DFT as linear transformation: 



The formula for DFT and IDFT is given by 
X(k) = Y J x(n)w k N n ,k = 0, 1 N-l (1) 



•(2) 



1 N-l 

t{n) = —Tx(k)w- N kn n = 0, 1, N-l 



Where W N =e- j27t/N (3) 

Let us define N-point vector x N of signal sequence x(n) , n = 0, 1, 
frequency samples and N*n matrix W N as 



.N-l and N-point vector X n of 



Xn- 



x(0) " 




X(0) 




x(D 


^w - 


X(l) 


w N= 


x(n— 1) 




X(AT-l) 





1 w 1 


w 2 


w N - 1 


w 2 


w 4 

"ft 


W 2(N-l) 



1 w 



N-l T i/2(JV-1) Ti/CAf-lKJV-l) 



N 



w, 



N 



w, 



N 



With there definition N-point DFT may be expressed in matrix 

X N = W N x N (4) 

Where, W N = matrix of liner transformation. 

Equation (4) can be inverted by per multiplying both side by W N " . Thus we can obtained x N = w N _1 X N 

(5) 

Which is expression for IDFT 

x n =1/NWn x X n (6) 

Where, W N = complex conjugate of matrix W N 

Wn' = W n 7n (7) 

Q. Compute DFT of 4-point sequence x(n) = { 0, 1, 2, 3} 
Solution: - 

X{k) = Y J x{n)e- i2AnlN 
X{k) = Y j x{n)e- jnknl2 

3 

X(0) = £ x(n)e° = x(0) + x(l) + x(2) + x(3) = 6 

n-0 

3 

X(l) = Y J x{n)e- im ' 2 = x(0) + x(l)e- j *' 2 + x(2)e- j * + x(3)e- j3 *' 2 = 6 

X(2) = -2 
X(3) = -2.4 

By matrix method, 





1 


1 1 


1 








w 4 = 


1 -7 -1 

1 -1 1 


j 
-1 






1 J - 1 


-J 






"0" 




'X(0)~ 




x N = 


1 

2 
_3_ 


x 4 = 


X(l) 

X(2) 
X(3)_ 




X 4 =W 4 * x 4 






"1 1 1 


1 " 




"0" 




= 


1 -j -1 
1 -1 1 


j 
-1 


* 


1 

2 






; j - 1 


-j. 




3 






0+1+2+3 


" 


5 




O-j-2 + 3; 




-2 + 27 


= 


0-1+2-3 




-2 




04 


j~2-3j_ 




-2 


-2 


7_ 



Q. Compute ID FT for the frequency component X(k) = { 60, 0, -4, 0} 

Solution: - 

x N =( W N */N ) X N x 4 =( W 4 */4) X 4 



X 4 



60 



-4 




w 4 



1111 

1 -j -1 j 

1-11-1 

1 j -1 -j 



x 4 =( W 4 /4) X 4 

"11 1 1 

1 -j -1 j 
1-11-1 

1 j -1 -j 



= Va 



60 


-4 




60-4 
60 + 4 
60-4 
60 + 4 



"56" 




"14" 


64 




16 


56 




14 


64 




16 



Properties of DFT:- 

The notation used to denote N-point DFT pair x(n) and X(k) as x(n) <- DFT-> X(k) 
1) Periodicity: 

Ifx(n)^DFT^ X(k) 

Where, x(n+N) = x(n) for all n. 

X(k+N) = X(k) for all k. 



Proof: 
x(n+N) = l/N£X(*)e' 2 * ( " 



+N)IN 



= \I^Y j X{k)e i2AnlN e i2A 



K=0 



= x(n) 
X(k+N)= ^n{k)e i2 ^ k+N)n 



11 = 

A'-l 



X(k)e J e J 



ii=Q 



=X(k) 



2) Linearity: - 

If xi(n) <r DFT^ Xi(k) and 

x 2 (n) <- DFT-> X 2 (k) Then, for any real or complex valued constants ai and a 2 
x(n) = aixi(n) +a 2 x 2 (n) <r DFT^ X(k) = aiXi(k)+a 2 X 2 (k) 
Proof: - 



aixi(n) +a 2 x 2 (n) <- DFT^ Y j a l x l {n) + a 2 x 2 (n) e ^ 1Anl1 



ii=() 



x x {n)e +a 2 2^x 2 (n) e 



ai X!(k) +a 2 X 2 (k) 
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Circular Symmetries of sequence: 



N-point DFT of finite duration sequence x(n) of length L < N is equivalent to N-point DFT of 
periodic sequence x p (n) of period N which is obtained by periodically extending x(n) i.e 

x p (n)= Y.xin-IN) 

x p (n) = x(n — k)=J] x(n -k- IN) 

J=-oo 

Finite duration sequence, 
x(n) = x p (n) 0<=n<=N-l 
= otherwise 

Is related to original sequence by x(n) by circular shift. This relationship is shown graphically as 
follows. 



This relationship is shown in figure 



x(n) 



12 3 



">•> 



">n 



-4-3-2-1012 3456 7 



x p (n-2) 



-6 -5 -4 -3 -2 -1 



2 3 4 5 



->n 



x'(n) 

3 



12 3 



">n 



x(2)=3 



x(l)=2 




x(0)=0 



x(3)=4 



x'(2)=l 



x'(l)=4 




x'(0)=3 



x'(3)=2 



In general the circular shift of the sequence is represented as index modulo N. i.e x(n) = x((n-k)) N 

For example, 

K = 2, N = 4 

That implies, 

x'(0) = x((-2)) 4 

x'(l) = x((-l)) 4 

x'(2) = x((0)) 4 

x'(3) = x((l)) 4 

Hence x'(n) is shifted circularly by 2 units in time, where the counter clock wise direction is selected as 

the +ve direction. Thus we conclude that circular shift of N- point sequence is equivalent to linear shift 

of its period extension and vise versa. 

Symmetry properties of DFT:- 

Let us assume that N-point sequence x(n) and its DFT are both complex valued. Then the sequence can 
be respresented as 

x(n) = x R (n)+jxi(n) 0<=n<=N-l (1) 

X' (k) = X R (k)+jXi(k) 0<=k<=N- 1 (2) 

X(k)==Y J x(n)e- j2m/N k = 0, 1, N-l 

«=o 

N-\ 

= V (x R (n) + jxj (n))(cos 2nkn I ' N - j sin lnkn I N) 

N-l N-l 

= ^ (x R (n) cos liikn I N + x t (n) sin 2nkn I N) - j V (x R (n) sin 2nkn I N -x } (n) cos 2nkn I N) 
= X R (k)+jX I (k) 

N-l 

X R (k) = = V (x R (n) cos 27ikn I N + x, (n) sin 27ikn IN) (3) 

N-l 

Xi(k) = Y, (x R O) sin 2nkn /N-x, (n) cos 2nkn IN) (4) 

«=o 

Similary, 



lV-1 



x R (n)=l/N ^(X s (fe)cos2^n/A^-X 7 (fc)sin2^n/A^) (5) 

N-l 

xi(n)= l/N^(X R (k)sm2nkn/N + Xj(k)co&2nkn/N) (6) 



k=0 



1) Real valued signal :- 

If x(n) is real, 

X(N-k) = X*(k) = X(-k) 

Proof: - 

N-l 

X(N-k)=Y j x(n)e- j2 * {N - k)n,N 



n=0 
N-l 



-jlm 



n=() 

=X(-k) = X*(k) 

2) Real and even sequence: - 

x(n) = x(N-n) 
Than, Xi(0) = , DFT reduces to 



X(k) =Y J x{n) cos iTtkn I N 0<= k<=N- 1 
IDFT Reduces to , 

N-l 

x(n) = 1/N Y, x (k)cos2nkn/N 0<=n<=N-l 



1=0 



3) Real and odd sequence: - 

If x(n) = x(N-n) 0<=n<=N-l 
Than, X R (k) = 0, DFT reduces to , 

N-l 

X(k) = -iY J x(n)sin2xkn/N ,0<=k<=N-l 

«=o 
IDFT reduces to, 

AM 

x(n) = j 1/N ^ X (k) sin 2nkn IN <= n<= N- 1 

k=0 

4) Purely Imaginary sequence:- 

x(n) = jxi(n) 

N-l 

X R (k) = V Xj (n) sin Ijtkn I N 

«=o 

N-l 

Xi(k) = ^ x, (n) cos Ijtkn I N 

If xi(n) is odd, then Xi(k) = and hence X(k) is purely real. 
If xi(n) is even then, X R (k) = and hence X(k) is purely imaginary. 



Multiplication of two DFTs and circular convolution :- 

Support we have two finite duriaotn sequences of length N, xi(n), X2(n). Their respective DFTs are, 



AM 



X 1 (k)=^jc 1 (n)e- j2 *' /iV k = 0, 1, N-l (1) 

n=Q 
N-l 

X 2 (k)=^x 2 (n)e- j2;An ' N k = 0, 1, N-l (2) 



»=o 



If we multiply two DFTs together the result in DFT say X 3 (k) of a sequence X3(n) of length N. Let us 
determine the relationship between x 3 (n) and the sequence xi(n) and x 2 (n) . 
We have X 3 (k) = Xi(k)X 2 (k) 



A'-l 



x 3 (n)= 1/N Y J ^ i (k)e i2Anlm 

k=0 
N-l 

x 3 (n)= 1/N Y J X i( k ) x i( k ) eilAnl " 



/.=() 



= i/nE 



x x (n)e 



/t=0 


L«=o 


A'-l A'-l 

= i/Njx,wj;i 2 (i) 

n={) 1=0 


Now, 


A'-l 

Y j a k =N a = \ 

k=0 


\-a N 
a- a 


vh* rf ± „ - J 2ll{m - 


n-l)/N 



j> 2 (/)^ 2M,// 



.1=0 



-jlnknlN 



2>~' 



- jlA(m-n-l) I N 



k=0 



N-l 

^a k = N -I = m-n+ pN = {{m-n)) N p is integer 



k=0 



= , otherwise 



Hence, 



x 3 (m)= 1/N Y J x 1 (n)x 2 ((m-n)) N N m = 0,l, N-l 

n=0 
N-l 

x 3 (m)= ^x 1 (n)x 2 ((m-n)) A , 



n=0 



The expression has the form of convolution sum. The convolution sum involves index ((m-n))N is called 
circular convolution. Thus we conclude that multiplication of DFT of two sequences is equivalent to the 
circular convolution of two sequences in time domain. 

Q. perform the circular convolution of the following two sequences. 
Xl (n) = {2,1,2,1} 

x 2 (n) = {1,2, 3, 4} 



Solution, 



N-\ 



x 3 (m)= ^x l (n)x 2 ((m-n)) 1 



n=0 



N = 4, 



*3(m)= ^x 1 (n)x 2 ((m-n)) 4 

n=() 
3 

x 3 (0)= ^x 1 (n)x 2 ((-n)) 4 =14 



xi(l)=l 



xi(2)=2 



x 3 (l)= ^x 1 (n)x 2 ((l-n)) 4 =16 




±> 



x 2 (l)=2 



x 2 (0)=l 





x 3 (2)= £jc 1 (n)jc 2 ((2-n)) 4 =14 



«=o 




x 3 (3) = ^ x l (n)x 2 ((3 - n)) 4 = 16 



^ 




«=() 




^ 




x 1 (n)x 2 ((3-n)) 4 =16 



x 3 (n)={ 14,16,14,16} 



Solution by DFT and IDFT method: 



Xi = 



x 2 = 



1 1 1 

-7 -1 j 

-1 1 -1 

7 "I -j 

1 1 1 ' 

-j "I 7 

-1 1 -1 

7 "I -7 





2 




6 


* 


1 









2 




2 




1 








1 




"10 


2 




-2-2; 


3 




-2 


L 4 




_-2-2;_ 



x 3 = 



60 





4 







IDFT: 

x 3 = (W 4 */4 )* X 3 



1 1 

1 -j 

1 -1 

1 j 



1 


1 " 


"60 " 


-1 


j 





1 


-1 


-4 


-1 


-j_ 






Parseval's theorem:- 

For complex valued sequence x(n) and y(n) 
x(n) <- DFT -^ X(k) 
y(n) <- DFT -> Y(k) 
Then, 

^x(n)y*(n) = ^X(k)Y*(k) 

Proof: - 



A'-l 



V x(n)y * (n) = r^ (0) Circular cross correlation sequence. 



«=o 



N k=Q 



1 Af-1 

r^'(0 = — Y J X(k)Y*(k)e j2M/N 

N k=o 



r xy (0) = ±-^X(k)Y*(k) 

N k=Q 

In Special case, y(n) = x(n) 

A<-1 JV— 1 1 JV— 1 

Y J x(n)y*(n) = Y J \x(n)\ 2 =-Y J \X(k)\ 2 

n=Q «=() iv n-0 

Which expresses the energy is finite duration sequence x(n) in term of frequency component {X(k)} 
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Fast Fourier Transform (FFT):- 

N-l 

X(k) = Y J x(n)e j2 ' !k,N 

The complex multiplication in direct computation of DFT is N and by FFT complex multiplication in 
N/2 log2 N . When number of points is equal to 4, the complex multiplication in direction computation of 
DFT is 16 and for FFT its value is 4. Hence the increment factor is 4. 



Radix-2 FFT algorithm :- 

If x(n) be the discrete time sequence then its DFT is given by X(&) = Vje(n)w^" k = 0, 1,....N-1 



N-l 



«=() 



Divide N-point data sequence into two N/2 data sequence fi(n) and f2(n) corresponding to the even 
number and odd number samples of x(n) respectively. 

fi(n) = x(2n) 

f 2 (n) = x(2n+l) n = 0, 1, N/2-1 

Thus fi(n) and fi(n) are obtained by decimating x(n) by a factor of 2 and hence the resulting FFT 
algorithm is called decimation in time algorithm. 



X(k)= j>(n)<* + j>(n)w 



n even 
N/2-1 



n odd 



X(k)= Y J x(2m)w k N 2 " + ^x(2m + l) 



w 



k{2m+\) 
N 



m-0 



m— 



But 



w 2 = ^-j2, nl 2j =e -j2,Kn / 2 )=WN ^ 



NI2-X 



NI2-X 



X(k)= YfM^Nn+WN ZA^K^ 
m-0 m-0 

X(k) = F 1 (k)+W„ i F 2 W 

Where Fi(k) and F2(k) are N/2 point DFT of sequences fi(m) and f2(m) respectively. 



x(0)- 
x(2)- 
x(4)- 
x(6)- 



x(D- 
x(3)- 
x(5)- 
x(7)- 



N/2 point 
DFT 



Fi(0) 



N/2 point 
DFT 




*X(7) 



Having performed that DIT once, we can repeat the processor for each of sequence fi(n) and f 2 (n) 

Thus fi(n) would result in two N/4 point sequence. 
v 11 (n) = f 1 (2n) 

vi 2 (n) = fi(2n+l) , n = 0, 1, N/4-1 

And f 2 (n) would result. 
v 2 i(n) = f 2 (2n) 

v 22 (n) = f 2 (2n+l) , n = 0, 1, N/4-1 

Then 



Fi(k) = vn(k)+ w N/2 v n (k) 



x(0)- 

x(D- 

x(3)- 
x(4)- 

x(5)- 
x(6)- 

x(7)- 



X m (p) 



X m (q) 



X 



W ° 



W 4 

VV N 



"X^ 7 



W 4 



X 



w ° 

VV N 



W 4 



X 



W ° 



W 4 




w ° 



w 2 

VV N 



W 4 



W 6 




w ° 

VV N 



W 2 



W 4 



W 6 




Fig: 8-point DIT FFT algorithm 

X m+ i(p) X m (p> — 



x m+i(q) 



X m (q)—^ 



w N ° 




♦ X(0) 



X(l) 



♦ X(2) 



♦ X(3) 



♦X(4) 



W 5 



♦X(5) 



W 6 



♦ X(6) 



W 7 




X(7) 



X m+ i(p) 



x m+ i(q) 



X m+1 (p) = X m (p)+W N r X m (q) 



Fig: Basic butterfly computation in DIT FFT 

algorithm. 



r+N/2 



X m+ i(q) = X m (p)+W N r+1N,z X m (q) 
= X m (p)-W N r X m (q) 



[W 



N/2 



■1] 



# Compute 8 point DFT for the sequence x(n) = { 1,2, 3, 4,5 ,6 } . Using DITFFT algorithm or DIF 
FFT algorithm. 

x(0) -^ j? ^ j, ^ z* • X(0) 





w 2 

VV N 



W 3 

VV N 



_J^. 



w ° 



-> 



w 2 




-> 



w ° 

VV N 



-> 



W ° 

" N 



-> 



W ° 

" N 



-> 



w ° 

VV N 



-> 



w ° 

VV N 



>X(1) 



-»X(2) 



-»X(3) 
-»X(4) 
-*X(5) 
-»X(6) 
-»X(7) 
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Q. Compute 4-point DFT for the sequence x(n) = { 14, 16,14, 16} using FFT algorithm. 

28 28 

x(0)=14 

x(2)=14_> 
x(l)=16 



x(3)=16-> 




Q. Compute 4 point DFT for the sequence x(n) = { 14, 16, 14, 16 } using FFT algorithm. 

x(0) ^- -J, ^— — ^ . X(0)=60 

*X(1)=0 

* X(2)=-4 
« X(3)=0 




Chapter: - 3 

Z-transform:- 

The z-transform of a discrete-time signal x(n) is defined as the power series. 

X(z)= f>(n)z- (1) 

When z is complex variable. 
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-3 



The inverse procedure( i.e obtaining x(n) from X(z) ) is called inverse z-transform . z-tranasform of a 
singal x(n) is denoted by 

X(z) = Z{x(n)} (2) 

Where as the relationship between x(n) and X(z) is indicated by 
x(n)^z -»X(z) (3) 

ROC (Region of convergence): 

ROC of X(z) is the set of all values of z for which X(z) attins a finite values. 

# Determine z-transform of following finite duration signals. 
xi(n) = { 1,2,5,7,0,1} 

X(z)= f>(n)z-« 

n-— °° 

= l+2z _1 +5z" 2 +7z" 3 +z" 5 
ROC: entire z-plane except z = . 

(2) x 2 (n)={2,4,5,7,0, 1} 

X 2 (z)= X(z)= X^ 2 (")z"=2z + 4z 2 +5 + 7z- 1 +z 

n— —oo 

ROC: Entire z-plane. Except z = and z = °° 

(3) x 3 (n) = 5(n) 
X 3 (n)= f> 3 (n) Z -" 

n—— oo 

= [S(n) Z -"l_ 
= 1 
ROC: Entire z-plane. 

(4) x 4 (n) = 5(n-u) k > 

X(z)= j>>4<X>z~" =2z + 4z 

X.4(Z) = B=-oo 

[S(n-k)z-t- k 
ROC: Entire z-plane except z = 

(5) x 5 (n) = 5(n+k) k>0 
X 5 (n)= X(z)= Y J x 5 {n) z - n 

n—--=<> 

= [S(n-k) z -"] n __ k 

-k 

= z 

ROC: Entire z plane except z = 0. 
ROC relationship to casusality:- 



Let us express complex z-variable in polar form as. 

z = re je (1) 

Where r = Izl = angle z . Then X(z) can be expressed as, 

X(z)l z=reje = X(z)\ = £*(n)r-V* (2) 

In the ROC of X(z) . \X(z)\ < °° 
But, 

|X(z)|= Y;X{n)r- n e- jen 

<2|jc(n)r-| (3) 

Hence X(z) is finite of the sequence x(n)r~ n is absolutely summable. 



The problem of finding ROC for x(z) is equivalent to determining the range of values of r for which the 
sequence x(n) r" n is absolutely sum able. 



|X(z)|< 2kn)r""|+Jkn)r"" 



n-0 



|X(z)|<£k-n)r-|+£ 



«=o 



x(n) 



•(4) 



If x(z) converges in some region of complex plain both summation in equation (4) must be finite in that 
region. 

If the first sum in equation (4) converges their must exists values of r small enough such that the 
product sequence x(-n) r 11 , 1 <= n < oo is absolutely summable. Therefore ROC for the first sum 
consists of all points in a circle of some radius ri where ri < oo as illustrated in the figure. 

V JZ) 
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>nd . 



Now if the 2 term in equation (4) converges there must exists values of r large enough such that 
x(n) 



product sequence 



hence ROC for second sum in equation (4) consists of all points outside a 



circle of radius r >t2 as shown in fig. 



W 




Since the convergence of X(z) requires that both sums in equation (4) be finite it follows that ROC of 
X(z) is generally specified as the annual region in the z plane r 2 <r<ri which is common region where 
both sums are finite. Which is shown in figure. 



R c (z) 




If r2 >ri there is no common region of convergence for the two sums and hence X(z) does not exist. 




Numerical: 

# Determine z-transform of the signal. 
X(n) = a n u(n) = a n , n => 
=0 , n<0 



Solution: 



X(z)= £jc(n)z"" 

= f>V" 

n=—oo 



1 



l-«z 
ROC: Izl > la I 



— for \az < 1 



x(n) 
4\ 



-> 



# Determine z-transform of the signal. 
X(n) = -a n u(-n-l) = -a n , n<=-l 

=0 , n>0 
Solution: - 

X(z)= £jc(n)z"" 

=-£(«■' V 

= -^ d (a~ l zj =-[{a- l z) + {a- l zf+. 
for \a~ x z\ < 1 






\-a l z 



i m (z) 




\-az 



ROC: \z\<\a\ 



x(n) 
-1\ 



-3 -2 -1 



-> 



IJz) 




R c (z) 



# Determine z-transform of the signal. 
X(n)= a n u(n)+b n u(-n-l) 
Solution: 



X(z)= £ac(n)z-" 

-1 

b z +2^a z 

= ±(b- l z) l +±(az- l y 
i= 

st p< 

i,e Izl > lal 



1=1 



n=0 



The first power series converges if lb _1 zl < 1 i.e Izl < Ibl and second power series converges if laz" 1 ! < 1 



Case I: Ibl < la I 



i m (z) 




R c (z) 



Case II: Ibl > la I 



Uz) 



X(z) = 



■ + - 



\-bz~ 1 l-az 1 
b-a 



a + b- z-abz' 
ROC: lal<lzl<lbl 




R c (z) 
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(1) linearity : 

If,xi(n) ^z^Xi(z) 

, x 2 (n) ^ z ^ X 2 (z) 

Then x(n) = aixi(n)+ a 2 x 2 (n) = aiXi(z)+aiX2(z) 
Proof: 

X(z)= £jc(n)z-" 

= V (a^j (n) + a 2 x 2 (n))z~" 

n—-c<> 

= a x ^x 1 (n)z _ " +a 2 ^x 2 (n)z~ n 
= a!Xi(z)+aiX 2 (z) 



# Determine z-transform and ROC of the signal 

x[n] = [ 3(2 n )-4(3 n )] 
Proof: 
Let, xi(n) = 2 n u(n) 
x 2 (n) = 3 n u(n) 
Then, 

x(n) = 3xi(n) -4x 2 (n) 
Accordign to linearity property. 



X(z) = 3Xi(z)-4X 2 (z) 
X 1 (z) = Z{x 1 (n)} = 



1 



X 2 (z) = Z{x 2 (n)} 



1-2Z- 1 
1 



ROC: Izl > 2 
ROC: Izl > 3 



l-3z 

The intersection of ROC of Xi(z) and Xi(z) is Izl > 3. 
Now 

3 4 



X(z) 



1-2Z" 1 l-3z _1 



Izl > 3 



# Determine z-transform of the signals: 

(a) x(n) = (cosw n) u(n) 

(b) x(n) = (sinw n) u(n) 

Solution: 

x(n) = (cosw n) u(n) 

= l/2(e iw " n +e- iw " n )u(n) 

(\ ■ 1 "i 

U 2 J 

Using linearity property. 

X(z) = -Z{e iw ° n u(n)} + -Z{e- jw ° n u(n)} 



Z{e jw " n u(n)} = 



1 



1 



Z{e' JW °u(n)} = 



l-e jw ° n z~ l \-2z~ l 
1 1 



1- TWeJl —1 1 O -1 

-g J " z l-2z 



ROC: Izl > l 
ROC: Izl > l 



X(z)=i 



l 



+ - 



l 



l 



2 1 jW,,n -l O 1 _ JWrJl 

\-e J a z 2 1-r'z 



1 - z Cosw,, 



l-2z 'cosw + z 2 
1 



ROC: Izl > 1 



(b)x(n)= —(e JW ° n e- JW " n )u(n) 

2; 



X(z) = 



1- z l Sinw Q 



1 - 2z cos w n + z 



ROC: Izl > 1 



(2) Time Shifting: - 

If x(n) <r z -> X(z) 
Then , x(n-k) <- z -> z" k X(z) 

Proof: 

x(n-k)^z^ j^x(n-fc)z~" 



Put n-k = m 



^x(m) 



-(m+k) 



= z~ k £x(m)z- (m) 



z" k X(z) 



# By applying the time-shifting property determine z-transform of x 2 (n) and x 3 (n) form z-transform of 
xi(n) given as, 

x x (n) = { 1,2,5,7,0,1}, x 2 (n)={l,2, T , 7, 0, 1} 
x 3 (n)= {0,0,1,2,5,7,0,1} 
Solution: 

X 2 (n) = xi(n+2), x 3 (n) = x x (n-2) 
Now, 

X 2 (z) = z 2 Xi(n) 
Xi(z) = l+2z" 1 +5z" 1 +7z" 3 +z" 5 
X 2 (z) = z 2 +2z+5+7z _1 +z" 3 



ROC: Entire z-plane except z = and z = oo 
X 3 (z) = z^XjCz) 

= z" +2z" +5z~ +7z" +z~ 
Roc: entire z-plane except z = 0. 



# Determine the z-transform of the signal . 

fl 0<n<A^-l 

x(n) = < 

[0 otherwise 

Solution: 

N-\ 



X(z)=2l. z -" 



l + z-'+z-' + z-' + z- (N - 1} =N ifz=l 

_, if z*l 



-, -N 

\-z 



\-z 



Alternative method :- 

X(n) = u(n) - u(n-N) 
Using linearity and time shifting property. 
X(z) = Z{u(n)}-Z{u(n-N)} 

\ z 

\-z \-z 

1 — ~ N 
= Z — ROC: Izl > 1 

l-z- 1 
Next method 



u(n) 












c 


o c 


( 


c 
















— > 



u(n-N) 
1\ 



N-l N 

(3) Scaling in z-domain:- 
If x(n) <rz -* X(z) ROC, ri <lzl <r 2 
Then a _1 x(n) <r z^ X(a _1 z) ROC: lalri <lzl <lalr 2 
For any real constant a real or complex. 

Z[a n x(n)]= '^ j a"x{n)z " = ^ x(ri){a~ n z)" = X(a~ 1 z) 

ROC : ri<la _1 zl <r 2 

lalrj < Id < \a\r 2 



N 



-> 
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# Determine z-transform of the signals. 

(a) x(n) = a n (cosw n)u(n) 

(b) x(n) = a n (sinw n)u(n) 



Solution: 

x(n) = a n (cosw n)u(n) 
xi(n) = cosw n u(n) 

Xi(z) = Z[cosw nu(n)] 



1 - z cos w 
l-2z _1 cosw> + z _1 



ROC: Izl >1 



Then, 
X(n) = a n Xl (n) 

X(z) = Z[a" Xj(n)] = X^a^z) = 



1 - az cos w 



1 - 2az L cos w Q + a z l 



ROC: Izl >lal 



(b) x(n) = a n (sinw n)u(n) 



X(z) = 



az x sin w. 







1 - 2az ' cos w n + a z 



ROC: Izl >lal 



(4) Time reversal: 

If x(n) <rz -» X(z) ROC, t x <lzl <r 2 
Then x(-n) <r z^ X(z _1 ) ROC: l/r 2 <lzl <l/ri 



Proof: 



Z[{x(-n)}]= ^(-»)z" 

n--oo 

Put , 1 = -n 
= f>(Qz' 



/=- 



= £*(0(z -1 )- 



= XCz" 1 ) 

ROC of X(z _1 ) is n <lzl <r 2 or, l/r 2 <lzl <l/n 
Note that ROC of x(n) is the inverse of that for x(-n) 



# Determine the z-transform of the signal. 
x(n) = u(-n) 
Solution: 

u(n) <rz^ — l — Izl > 1 
1-z 
Using time reversal property. 

u(-n) <- z -> -^— Izl < 1 
1-z 



(5) Differentiation in z- transform 

If x(n) <r z -» X(z) 

T , M , v dX(z) 
Then n x(n) <r z -r -z 

Proof :- 



X(z)= £^(n)z-" 

n—— °° 

rentiation, 

= ^*(n)(-n)z~" 



By differentiation, 
<ff(z) 



dz 

= -z~ 1 ^{«x(«)z~" 

= -z~ 1 Z{nx(n) } 

Z{nx(n) } = -z Note that both transform have same ROC. 

# Determine z-trnasform of the signal. 

x(n) = na n u(n) 
Proof: 

xi(n) = a n u(n) 



xi(n) 



a — az 



ROC; Izl >lal 



x(n) = n x^n) 
X(z) = Z{n Xl (n)} 



dX x {z) 
dz 



■1 



= - z 



(l-az- 1 ) 2 



az 



-2 



az 



(l-az' 1 ) 2 



ROC: lzl>lal 



# Determine the signal x(n) whose z-transform is given by 
X(z) = logCl+az" 1 ), lzl>lal 



Solution: 

dX(z) 1 



dz 
dX(z) 



a + az' 1 ) 



az 



az 



dz (1 + oz ) 
1 



■az 



}-{-a)z- 
Taking inverse z-transform. 



n-l 



nx(n) = a(-af u(n-l) 
x(n) = (-l) nl . a n /n u(n-l) 



Date: 2066/08/26 



If xi(n) <-z -> Xi(z) 
x 2 (n) <r z^> X 2 (z) 

Then 

X(n) = xi(n).x 2 (n) <rz -> X(z) = Xi(z) X 2 (z) 
The ROC of X(z) is at least the intersection of the for Xi(z) and X 2 (z). 

Proof: 

The convolution ofxi(n) and x 2 (n) is defined as 



x(ri) = y]x l (k)x 2 (n — k) 
z-transform of x(n) as, 
X(n)= f>(n)z-' 



X{n)= £ Y 4 x x (k)x 2 {n-k)z- n 

n--oa\_k--ao 

X(n)= ^ Xl (k)^x 2 (n-k)z- n 
X(n)= ^x I (/r)z""x 2 (z) 



k=- 



X{n) = X 2 {z)Y J x l {k) Z - n 



k=- 



= X 2 (z)X l (z) 

# Compute convolution x(n) of the signals 
X(n) = { 1,-2,1} 

\1 < n < 5 
x 2 (ri) = \ 

otherwise 



Inverse z-transform. 



X(z)= £x(*)z"* (i) 



k=- 



n-l 



Multplyign both sides of (i) by z " and integrate both sides over a closed contour within the ROC of 
X(z) which enclosed the origin . 

Now, 



§X(z)z"- l dz = <f ^xikW^dz (2) 



Ck=- 



Where C denotes closed counter in the ROC of X(z) . 

We can interchanged the order of integration and summation on right hand side of (2) . 



| X(z)z n - 1 dz = ^x(k)$z n - l - k dz (3) 



k=-<* 



From Cauchy integral theorem , 



Figure: 



— <pz" 


dz = 


Jl k = n 
[0 k*n 


Then, 







§X(z)z"-'dz = 27gx(n) 

x(n) = ix{z)z"- l dz (4) 

2 ^'c 



7) Multiplication of two sequences: 



^ 



v l dv 



If xi(n) <-z -» Xi(z) 
x 2 (n) ^ z^ X 2 (z) 
Then 

x(n) <- z^ xi(z)x 2 (z) <- z^X(z) = ix x (u)X. 

Where C is sthe closed counter that encloses the origin and lies within the ROC of common to both 
Xi(v)andX 2 (l/v) 

X(z)= f>(tt)z- 



\vj 



X(z)= ^X!(n)x 2 (n)z " 

n—— °° 

Where, 

x 1 (n) = —lx i (v)v"- 1 dv 
2 W c 



x 1 (n) = -^ T |x i (v)X 2 * 

2 ^c 



/ »\ 



v dv 



\ v J 



X(z) 



= ^rl^i(v) 



2^ J 



X X 2(") 



^V» 



vv; 



v rfv 



8) Parse val's theorem :- 

If xi(n) and x 2 (n) are complex valued sequence, then, 



1 



Y J X l (n)x 2 *(n) = — \x i (v)X t 2 (\lv')v- l dv 
2xj J c 



X(n) = c n for all n. 

Where X(z) is rational, the expansion can be performed by long division. 



# Determine inverse z-transform of 

X(z)= l 



l-0.5z-'+0.5z- 2 



When 

(a) ROC: Izl > 1 

(b) ROC : Izl <0.5 



Solution: 
(a) ROC: Izl >1 , x(n) is causal signal. 



l-0.5z" 1 4fl.5z" 2 ) 1 (l-O.Sz^-OMz" 

-1 -2 

l-0.5z +0.5z 

~1 Zl 

0.5z -0.5z 

-1 -2 -3 

0.5z -0.25z +0.25z 

71 7T~ 

-0.25z -0.5z 

-2 -3 -4 

-0.25z +0.5z -0.125z 

3 31 

0.375z +0.125z 

-3 -4 

-0.375z +0.18175z +0.18175z 

3J 3 

-0.0625z +1.8175z 



X(z) = {l+O.Sz^-O^Sz^-O.STSz" 3 } 

x(n)={ 1,0.5,-0.25,-3.75, } 

(b) ROC : Izl <0.5 x(n) is anticausal signal. 



-2 -1 

0.5z -0.5z +1) 1 



2 3 4 
(2 z + 2z - 2z 



l-z+2z 




2 
z - 2z 

2 3 
z - z + 2z 




2 i 
z - 2z 

2 3 4 
- z + z - 2z 




3 4 

-3z 3+ 2z 4 

-3z +3z - 


5 
6z 


4 


i 



-z + 6z 



X(z) = 2z 2 +2z i -2z 4 -6z 5 

= 6z 5 -2z 4 

.-.*(*) = { -6,-2,2,0,0} 



6z 5 -2z 4 +2z 3 +2z 2 +0 + 



Initial value theorem: 

X(0) = lim z -> oo X(z) 

Final value theorem: - 

Lim n— > oo x(n) = lim z— >1 (z-1) X(z) 



# The impulse response of relaxed LTI system is h(n) = a"u(n) \a\<l. Determine the step response 

of the system where n tends to infinity. 
Solution: - 

y(n) = x(n)*h(n) 
x(n) = u(n) 
h(n) = a n u(n) 
Y(z) = X(z)H(z) 

r i V i "i 



X(z) 



1-z" 



l-az~ 



Now, 
Lim n —> go y(n) = lim z —>l (z-l)Y(z) 

2 

= lim z tends to 1 (z-1) 

(z-D(z-a) 

1 



1-or 



System function (transfer function) of LTI System: - 

The o/p of LTI system to an input sequence x(n) can be obtained by comparing the convolution of x(n) 
with unit sample of the system, i.e y(n) = x(n)*h(n) 
We can express this relationship in z-domain as 
Y(z) = X(z)H(z) 

When, Y(z) = z-transform of y(n) 
X(z) = z-transform of x(n) 
H(z) = z-transform of h(n) 

XT XX, N 7(Z) 

Now, H(z) = 

X(z) 

H(z)= f>(n)z"" 

n— -co 

H(z) represent the z-domain characteristics of the system where as h(n) is corresponding time domain 
characteristics of the system. The transform H(z) is called the system function or transfer function. 

General form of liner constant coefficient difference equation: - 

The system is described by liner constant coefficient difference form. 

N M 

y(n) = -£ a k y(n-k) + Y J b k x(n-k) (1) 

*=1 k=0 

Taking z transform on both sides, 



Y(z) = -Y j a kZ - k Y{z) + Y.b kZ - k X{z) 





k=\ 


k=() 


Y(z) 
H(z 


r n \ 


M 

= Yb k z- k X{z) 

k=Q 
M 

lhz- k 

=0 


X(z) , , 


N 

V -* 



•(2) 

^ a kZ'" 

Therefore a LTI system described by constant coefficient difference equation has a rational system 
function from this general form we obtained two important special forms: 

M 

1. if ak = , for 1 <= k<= N . Equation (2) reduces to H(z) = "5]b k z~ k (all zero system) such a 

*=o 
system has finite duration impulse response and it is called FIR system or moving average 
(MA) system. 



2. On the other hand if bk = , for 1 <= k<= n the system function reduces to 



H(z) = jt^- (4) (all pole system) Due the presence of poles the impulse 

-k 

k=l 



k=l 

response of such a system is infinite in duration and hence it is IIR system. 



# Determine the system function and unit sample response of the system describe by the difference 

equation y(n) = Vi y(n-l) +2x(n). 

Solution: 

Taking z-transform on both sides, 

Y(z) = l/2Y(z) + 2X(z) 

(l-l/2z- 1 )Y(z)=2X(z) 

HW =2M=_?-_ 

X(z) \-\l2z l 
By inversion, 
h(n) = 2 (l/2) n u(n) 
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Response of pole zero system with Non-zero initial condition:- 

The difference equation, 

N M 

y(n) = -£ a k y(n-k) + ^b k x(n-k) 0) 

k=\ k=Q 

Suppose that the signal x(n) is applied to the pole zero system at n = . Thus the signal x(n) is assume to 
be a causal. The effects of all previous input signal to the system are reflected in the initial condition y(- 

1), y(-2) y(-N). Since the input x(n) is causal and since we are interested in terminating the o/p 

signal y(n) for n > we can use one sided z-transform which allows us to deal with initial condition. 

Now, 



Y + {z) = -Y j a kZ - x 



k=\ 



k=0 



Since x(n) is causal, We can set X + (z) = X(z) 

M N k 

Yj>kZ k Yj a ^ k Yjy(- n ^ z ~ n 

Y + (z) = ^ .X(z)-^ f 5 

k-\ k-\ 

N (z) 
= H(z)X(z) + -^ (3) 

Mz) 

Where, 

N (z) = -f j a kZ - kJ £y(n)z" (4) 

k=\ n=\ 



III 

+ Yb kZ - k X + (z) (2) 



From (3) the o/p of the system can be sub divided into two parts. The 1 st part is zero state response of the 

system defined in z-domain as Y zs {z) = H{z)X{z) (5). 

The second component corresponding to o/p resulting form initial condition. This o/p is zero input 

response of the system which is defined in z-domain as Y .■ (z) = — (6) 

Mz) 

Since the total response is the sum of the two o/p component which can be expressed in time domain by 
determining inverse z-transform of Y zs (z) = Y zi (z) . Seperately and adding the result. 
y(n) = y zs (n) + y zi (n) (7) 

# Determine the unit step response of the system described by the difference equation y(n) = 0.9y(n-l) - 
0.81 y(n-l) + x(n). Under the following initial condition. 

(a) y(-l) = y(-2) =0. 

(b)y(-l) = y(-2) = l. 

Taking one-side z-trnasfrom. 

Y + (z) = 0.9z- l [Y + (z) + y(-l)z\-0.8z- 2 [Y + (z) + y(-l)z+y(-2)z- 2 ]+X + (z) 
Fory(-l) = y(-2) = 

Y + (z) = 0.9z- 1 Y + (z)-0.8l Z - 2 Y + (z) + X + (z) 

Y + (z) = 1 

x + (z) l-0.9z '-0.81z 2 

Y:(z)_ z 2 z 2 



z (z z -0.9z + 0.8z)(z-l) (z-0.9e J *' J )(z-l) 



Z K-\ 1 2 



+ -^^ + " 



(z 2 -0.9e j ' rn )(z-0.9e- j ' rl3 )(z-l) (z-0.9e^ /3 ) (z-0.9e^ /3 ) z-1 



*!= 1.098c 



-75.53° 

*.= 1.098c'' 5 - 53 ' , k 2 = 1.1 



+ 1.098g' j553 1.098g A53 1.1 

wCz) "a-o.9^ / v i ) + a-o.9e-^ / v i ) + i-z- 1 

By inversion, 

y zs (n) = [l.098e" A53 (0.9e 7Vr/3 )" +1.098e j553 (0.9e"^ /3 )" +l.lj+ y(n) 
= [l.l + 1.098(0.9)" cos(n7T/3-5.53)\i(n) 

(b) y(-l) =y(-2) =1. 

Y Z M)= "^ 

A(z) 

N k 

N {z) = ^ k z- k ^y{n)z- n 

k=\ n=\ 



N (z) = -f j a k z- k Y j y(-n)z- n 

1 n=\ 

= [a iZ - l y{-\)z + a 2 z- 2 {y{-\)z+y{-2)z 2 )] 

N (z) = -[-0.9x1 + 0.8 lz _1 + 0.8 lxl] = 0.09-0.81Z -1 



Y,(z) = 



0.09-0.81Z _ 
l-0.9z _1 +0.8iz~ 2 



Y(z) = Y zs (z) + Y zi (z) 



1 .099 0.568 + y'0.445 0.568 - y'0.445 

■ + —7-, 7- + ■ 



\-z 



1-0.9^'V 1 



l-0.9e*'V 



y{n) = [l .099 + 1 .44(0.9)" cos(n^/ 3 + 38)]w(n) 
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# Determine well known fibanacci sequence of integer numbers is obtained by computing each term as 
the sum of two previous ones, the first few terms of the sequences are 



1,1,2,3,5,8 

Determine a close form expression for the nth term of Fibonacci series. 

Solution: - 

Let y(n) be the nth term. Then, 

Y(n) = y(n-l) + y(n-2) (i) 

With initial condition , 
y(0) = y(-l)+y(-2)=l 

y(D = y(0)+y(-D=i 

y(-l)= 1-1 = 0, y(-2)= 1 . 

Taking one sided z-tranzsform on both sides of (i). 

Y + (z) = z- l Y + (z) + y(-l) + z- 2 Y + (z) + Y(z) + yi-Vz' 1 

Y + (z) = z- l Y + (z) + z- 2 Y + (z) + l 




Taking inverse z-transformation. 



y(n) = 



1 + V5 

2V5 



1+VsV fi-VsYi-VsY 

2V5 



w(n) 



# Determine step response of the system y(n) = a y(n-l) + x(n) -1 < a < 1 

When initial condition is y(-l) = 1. 
Solution: 

Taking one sided z-transform, 

Y + (z) = a[z- 1 Y + (z) + y(-l)\+X + (z) 

X + (z) = z[u(n)] = — i -^ 
\-z 

Or, Y + (z) = az^Y + (z) + a + - 



\-z 



Or, Y + (z)(l-az 1 ) = a + - 



1 



Y + (z) = 



Y + (z) = - 



a 



■ + - 



\-z 

1 



l-az 1 (l-z^a-az' 1 ) 



a 



- + 



l-az 

Taking inversion, 

y(n) = aa"u(n).l 



\\-aj 



\-z 



f 1 A 



\\-ccj 



l-az 



\l-a y 



u(n)- 



\-a 



a"u(n) 



a n+l (l-a) + l-a 



y(n) = 



\-a 

l-a n+2 ~ 
\-a 



u(n) 



# Determine the response of the system y(n) = 5/6. y(n-l) - 1/6 y(n-2) + x(n) to the input signal 

x{n) = 8{n) 8{n - 1) 

Solution: 

Taking z-transform in both sides. 

Y(z) = ^z- 1 Y(z)-^z- 2 Y(z)+X(z) 
6 6 

X(z) = l--z l 



Y(z) = -z- l --z- 2 Y(z) + 



f 1 A 
1 — z 

3 



1 J -1 , - 
1 — z +-z 

6 6 



Y(z) = 



1 — z 



1 — z 



Y(z) = 



V 



J 



Y(z) = 



( 5 \ \ 

I 6 6 J 

1 



i4z- 



Taking inverse z-transform, 



y(n) = \ 



u(n) 



Causality and stability: - 

A causal LTF system is one where unit sample response u(n) satisfies the condition h(n) = 0, for n 
< 0. We have also shown that ROC of z-transform of a causal sequence is exterior of a circle. 

A necessary and sufficient condition for a LTI system to be BIBO stable is V \h(n)\ < °° . In turn this 

n=-oo 

condition implies that H(z) must contain the unit circle within it's ROC. 
Indeed, Since 

H(z)= f>(n)z" n 



\H(z)\< 



2>(*)z~" 



When evaluating on the unit circle i.e \z\ = 1 



H(z)< 



Y^Kn) 



Hence if the system is BIBO, the unit circle is contained in the ROC of H(z). 



# A LTI system is characterized by the system function. 



H(z) = 



3-4z~ 



l-3.5z^+1.5z" 2 



Specify the ROC ofH(z) and determine h(n)for the following condition. 

(a) The system is stable. 

(b) The system is causal. 

(c) The system is anticausal. 



Solution: - 
H(z) = 



1 — z 

2 



- + - 



l-Sz- 1 



The system has poles at z = Vi and z = 3. 



(a) Since the system is stable, its ROC must include the unit circle and hence it is Vi < Izl <3. 
Consequently h(n) is non-causal. 



h(n) = 



r j\" 



u(ri) - 2.(3)" u(-n - 1) ( The system is unstable). 



v^y 



(b) Since the system is causal , its ROC is Izl >3. In this case h(n) = 

system is unstable). 

(c) Since the system is anticausal it's ROC is Izl <l/2 . In this case, 



r i\" 



\^j 



u(n) + 2(3)" u(n) (The 



h(n) = 



u(—n - 1) + 2.(3)" u(-n - 1) . (The system is unstable). 
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Transient and steady state response: 



# Determine the transient and steady state response of the system characterized by the difference 



equation y(n) = 0.5y(n-l) +x(n). When the input signal is x(n) = lOcos 
initially at a rest (i.e it is relaxed ). 



nn 



u(n) . The system is 



Solution: 

Taking z-transform, 

Y(z) = 0.5z- l Y(z) + X(z) 

Y(z) = 1 

X(z) 1-0.5Z 1 

Z(cosw nu(n)) 



1 - z cos w 



1 - 2z cos w n + z 



X(z) = 10 



1 -i n 
\ — z cos — 



. - _! K _ 2 

l — 2z cos — \- z 



10 



V2"' 



y 



Fv? 



c j +z- 2 ) 



y(z) 



(l-0.5 Z - 1 )' JU pV2? r +7 r j 
1 



.10 
10(1 



V2" 



z" 1 ) 



(\-0.5z- l )(l-e 1 « l \- l )(\-e-""\- 1 ) 



-jit I A -In 



63 



6.78e- A870 6.78^' 28 ' 7 

+ r-; r + 



1 at -1 ' i i/r/4 -1 

1-0. 5z l-e J z 



i -;>/4 -1 

1-e z 



The natural or transient response is 

y nr (n) = 6.3(0.5)"w(n) and forced or steady response is 

y fr (n) = [6JSe mJ e jn ' r ' 4 +6JSe mJ e- jn '" 4 Un) 



= 13.56cos 
Pole-zero diagram: 
H(z) 



n — -28.7 

v 4 



u(n) 



l-z-'-2z- 2 



l-1.75z _1 +1.4z -0.375z~ 

3 2 

z -z -2z 
z 3 -1.75z 2 +1.25z-0.375 



z(z-2)(z + l) 



(z - 0.75)(z - 0.5 - ;0.5)(z - 0.5 + y0.5) 



<r 



-©- 



-l 



A 



J0.5 



I m (z) 



-X 



-$h 



-X- 



0.75 1 



2 R e(z) 



J0.5 X 



Fig: pole-zero diagram 



Zeroes at z = 0, 2, -1 
Poles at z = 0.75, 0.5±y0.5 



# A filter is characterized by following poles and zeroes on z-plane. 





Zeroes at 


Poles at 




Radius 


Angle 


Radius 


Angle 




0.4 


7irad 


0.5 


Orad 




0.9 


1.0376 


0.892 


2.5158 rad 


Shows z-plane plot and plot ma| 


mitude response (not to sea 


e). 


Solution: - 


Zeroes are, 


z = 0.4e J> =-0.4 


z = 0.9e i ' L0376 =0.45 + jO.77 


Poles are: 


z = 0.5e j0 =0.5 











z = 0.892e ZM35 = -0.728 + y'0.522 



<r 



A 



-0.722 



J0.77 



"©- 



-0.4 



Im(z) 



— 

jO.522 I 
I 
I 



0.45 



-* > 



0.5 



Fig: z-plane plot 



R e (z) 



H(z) = (z + 0.4)(z-0.45-j0.77) 

(z-0.5)(z + 0.77-y0.522) 

Put, z = e jwT , 

TT , , (e jwT + 0.4)(« jwT - 0.45 - 70.77) 
H(z) = 



JwT 



(e m - 0.5)(e Jwl + 0.77 - jO.522) 



H(e JW ' ) 



WT=Q 



\H(e ]W1 ) 



WT=x/4=0.319 



\H(e jm ) 



WT=3nl 4=3.438 
jwT - 



\H(e Jwl ) 



WT=n=l.U 



4 



HCei™ 1 ) 



0.478 




-> 



Ji/4 jt/2 331/4 71 



Fig: magnitude response 



coT 



Notch filter: 



It is a filter that contain one or more deep notches or ideally perfect nulls in its frequency response 
characteristics. Notch filter are useful in many application where specific frequency component must be 
illuminates for example, instrumentation in recording system requires that power line frequency of 50 
Hz and its harmonics be illuminated. 




COT 



To create a null in frequency response of filter at frequency wo we simply introduce a pair of complex 
conjugate zero on the unit circle, z, 2 = e 



:jw (l 



System function is given by , 

H(z) = b o a-e jw "z' 1 )a-e jWo z' l )(l-e jw "z' 1 ) 
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Response to complex exponential signal: Frequency response function:- 
The input-output relationship for LTI system is u(t) = V h(k)x(n -k) (i) 



&=-<* 



To develop a frequency domain characteristics of the system. 
Let us excite the system with complex exponentials. 



x(n) = Ae 



- A»J wn 



■ °° < n < °° 



•(2) 



Where, A is amplitude and w is orbitrary frequency confined to the frequency interval [-;r,;r] 
Now, 



y (n)= Y j h(k)[Ae Mn - k> ] 
^h(k)Ae- jwk 



k=-<* 



k=-<* 



, jwn 



•0) 



H(w)= ^h(k)e- jwk (4) 

The function H(w) exists if the system is BIBO stable, that is 
f>(/0|<~ (5) 



The response of the system to complex exponential is given by , 

y(n) = AH(w)e jwn (6) 

We note that the response is also in the form of complex exponential with the same frequency as the 
input but altered by multiplicative factor H(w). The multiplicative factor is called eigen value of the 
system in this case a complex exponential of the form (2) is an eigen function of LTI system and H(w) 
evaluated at frequency of input signal is corresponding eign value. 



# Determine the o/p sequence of the system with impulse response h(n) = 

complex exponential sequence x(n) = Ae iml2 u(n) - °° < n < °° . 

Solution: - 

H(z)= l 



1 



u(n) When the input is 



H(w) = 



i-V 

2 
1 

\--e~ jw 
2 



At w = 



n 



H 



'^ 



1 



-J26.6 



\^J , 1 -i\ V5 
1 e l 

2 
And therefore the o/p is 

\ . a 

e 2 



y(n) = A 



V?' 



-J26.6 



J 



2 K^-16.6) 

y(n) = ^Ae 2 



■ oo < n < °° 



# Determine the response of the system with impulse response h(n) = — u{n) to the input signal 

v2y 



71 

x(ri) = 10-5sin — n + 2cos7m 

2 



— °° < n < °° 



Solution: - 






H(z) = 




1 


1- 


1 i 

- — z 
2 


H(z)-- 




1 




--e~ jw 

2 




1- 



The first term in a input signal is a fixed signal component corresponding to w = 0. Thus, 



H(0) = 



= 2 



1 



The second term in x(n) has frequency — . Thus , 



H 



71 



1 



1 n 

1 e z 

2 



7F 



-j 26.6 



Finally the third term in x(n) has a frequency w = n . Thus H(;r) = 

Hence the response of the system to x(n) is, 

, , ™ 10 . f/r „, ^ 40 
y(w) = 20 = sin — n-26.6 H cosmi, —°° <n<°° 

S U J 3 



# A LTI system is described by following difference equation y(n) = ay(n - 1) + bx(n) 

0<a <1. 

(a) Determine the magnitude and phase of frequency response H(w) of the system. 

(b) Choose the parameter 'b' so that the maximum value of l/f (w)| and phase of H(w) for a= 0.9. 



it 



(c) Determine the o/p of the system to signal for x(n) = 5 + 1 2 sin — n - 20 cos 



mi-\ — 
v 4 
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Chapter:- 5 



Discrete filter structure :- 














Cascade form structure :- 












H(z) = Ut l H k (z) 












Where k is integer part of (N+l)/2 












„ //N b ko +b kl z~ 1 +b k2 z~ 1 
H k (k) = -* *!_ ^_ 

\ + a kl z +a k2 z 






x 3 (n) 
Y2( n ) 


x k (n) 








Hi(z) 


x 2 (n) 


H 2 (z) 


H k (z) 






Yl(n) 


yd 



x k (n) 




y k (n)=x k+1 (n) 



Parallel form structure: 



H(z 






H 



P k Z 




# Determine the cascade parallel realization for the system described by the system function 



10 



i-V 



H(z) = 



1 — z 



(l-2zi 



3 -i 



1 — z 

4 



1 i 



1 — z 
8 



1 + 



V 2 + i 2y 



'2 2 



Cascade form:- 



1 — z 



HAz 



1 -^ 



1--Z- 1 
3 

1 — z + — z 
8 32 



1 — z" 



(l + 2zi 



H 2 (z) = 




1 1 



1 -• i *- -2. 

\-z +-z 

2 

H(z) = lOH l (z)H 2 (z) 
The cascade relationship is 

x(n) 

> 1 




Parallel form: 



H(z) = 



10 



z — 



V 



z — 



{z + 2)z 



f 3V O 

z — z — 
I 4A 8 



2 + J 2 



2 J 2y 



#(z) 



10 



z — 



■J 



z — 



V 



fc + 2) 



•y 



l 



l 



v2 +7 2 y 



2" J 2 



- + - 



■ + - 



3 1 

z — z — z- 

4 8 



V 2 + J 2y 



+ - 



1 .1 



2 J 2j 



^=2.93, /t 2 =-17.68 , fc 3 = 12.25 -714.57, 



fc* =12.25 + 714.57 



2.93 -17.68 12.25+ /14.57 12.25+ /14.57 

H(z) = + : + ; : + ■ 



, 3 _, 1 _, 

1 — z 1 — z 1- 



1 .1 

2 + ; 2 



12 2 




Yk(n) 



Cascade form structure :- 



H(z) = n*H k (z) 



ko+ b ^ 



H k (z) = b kn +b uZ - l +b k2 z- 2 



k = l,2 



x(n) = x^n) 



Hi(z) 



H 2 (z) 



| H k (z) | > 

y(n) 



Xk(n) 



b t 



b,,2 



^o — <5 



V 

- 1 — > 



Yk(n) 



Lattice form structure: 

Fir filter with system function 
# m (z) = A m (z) 



m = 0,1 m — \ 



A m (z) = l + f j a m (k) Z - k m>\ 
k=\ 

A (m) = l 

h m (z) = l h m (k) = a m (k), k = l,2. 

m > degree of polynomial. 



.117 



y(n) = x(n) + £ a m (k)x(n - k) 



k=\ 



x(n) 



\/l 



Orn(l) 



a m (2) 



3k 3k 

k5 — <5 




V 



-> 



y(n) 



For m = 1 . 

y(n) = x(ri) + a x (l)x(n - 1) 

1 st order lattice filter. 



Figure: Direct form realization. 




x(n) 



go(n-D 



gl(n) 



f Q (n) = x(n) 

g Q (n) = x(n) 

fi(n) = f (n) + k l g Q (n-l) 

y(ri) = x(ri) + k x x(n - 1) 

g 1 (n) = k 1 f Q (n) + g Q (n-l) 
= k l x(n) + x(n-i) 

k x = reflection coefficient 




x(n) 



y(n) = x(n) + a 2 (l)x(n -l) + a 2 (2)x(n - 2) 
fi(n) = f x {n) + k 2 g x {n-l) 
= x(m + k x (1 + k 2 )x(n - 1) + k 2 x(n - 1) 
= a 2 (2) = k 2 , a 2 (l) = k l (l + k 2 ) 



Conversion of lattice coefficient to direct form filter coefficients :- 

Direct form filter coefficient {a m (k)} can be obtained form lattice coefficients {k i } using the relations. 

A Q (z) = B Q (z) = l 

A m (z) = Am-\z) + k mZ - x B m _ l (z) m = l,2, m-\ 

BJz) = z- n A m (z- 1 ) m = l,2, m-\ 

B m (z) -> Reverse polynomial of A m (z). 

# A given 3 -stage lattice filter with coefficient k, = — , k-, = — , k, = — . Determine FIR filter 

1 4 2 3 

coefficient for direct form structure. 

A l (z) = A (z) + k l z- l B (z) 

=i + V 

4 
a l (0) = 1, a x (1) = — , Corresponding to single stage lattice filter. 



1 

= — +z 



-i 



/3 l (z) = z- 1 A l ( z - 1 ) = z- 1 \l+ ] -z 

A 2 (z) = A l (z) + k 2 z- 1 j3 1 (z) 

= l + -z + — z 
8 2 

3 1 

a 2 (0) = 1, a 2 (1) = — , a 2 (2) = — Corresponding to 2 nd order lattice form. 

8 2 

A (z)= i + | z ' 1+z ' a 

A 3 (z) = A 2 (z) + k 3 z- 1 Mz) 

=1+ — z +-z +-z 
24 8 3 

13 5 1 

«r 3 (0) = l, a 3 Q) = —, a 3 (z) = - , a 3 (3) = - 

24 8 3 
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Conversion of direct form FIR filter coefficient to lattice coefficient: 

A m M)= AJZ \ +k ; 2 BJZ) m=m-l, -1 

K=oc m =a(m) 



# Determine lattice Coefficient Corresponding to FIR filter with system function. 



13 



1 



13 



H{z) = A i {z) = l + — z~ l +-z~ 1 +-z~ i + — z~ i + z 



24 



24 



Solution: 



k l =a 1 (3) = -, 



8 3 

B 3(z) = -+^Z^-r^ 7 Z^+Z 



1 . 5 _, 13 _ 2 , _ 3 



3 8 



24 



Mz) = 



A 3 (Z)-k 3 B,(z) .3 , 1 2 

— f-^ = l + -Z + — Z 

l-kj 8 2 



z -Z o 



A(z) 



A 2 (z)- J fi: 2 5 2 (z)_ 1 , 1 



k. = a, (1) = 



P^ = l + -r' 

1 



Lattice and lattice ladder structure for IIR system:- 

Let us begin with all-pole system with system function. 

H(Z0 = r— * = —^— CD 



l + ^a m (k)z 



-k 



Am(z) 



Difference equation for this system is 

N 

y(k) = ^a m (k)y(n-k) + x(k) 

WhenN= 1. 



(2) 




>y(n) 



x(n) = fKn) 

f (n) = f,(n)-k lgo (n-l) 
gi(n) = kif (n)+g o (n-l) 
y(n) = fo(n) 

= fi(n)+kig (n-l) 
y(n) = x(n)-kiy(n-l) 
yi(n) = kiy(l)+y(n-l) 



x(n) 



fo(n) 



f 2 (n) 



> 



g2( n ) 



gl(n-l) 



< 



(n) 





fo(n) 




y(n) 


>( x 


/ Vi 










Z" 1 






go(n-l) 




go( n ) 



Fig: Two stage lattice structure 



f (n) = x(n) 

fi(n)=f 2 (n)-k 2 gi(n-l) 

g 2 (n) = k 2 fi(n)-kig 2 (n-l) 

f (n) = f 1 (n)-k 1 go(n-l) 

gi(n) = kif (n)+g (n-l) 

y(n) = f (n) = g (n) 

= f 2 (n) - k 2 gi(n-l) - kig (n-l) 

= - ki (l+k 2 ) y(n-l) - k 2 y(n-2) + x(n) 



in 



-k 



H{z) 



k=0 



k=0 



CM 

KM) 



Ladder Co-efficient, 
v m =C m (m) m =0, 1,2, 3, 

C m - l (z) = c n (Z)-v m /B m (z) 

x(n) /^ f Q (n) 



.M 




# Obtain the lattice ladder structure. 

TT/ , l + 3z~' +3z~ 2 + z~ 3 

H(z) = i o T 

l-0.9z l +0.64z~ 2 -0.576z~ 3 

= C 3 (z) 
A 3 (z) 
Solution: 

Lattice coefficient:- 



A 3 (z) = l-0.3z -1 + 0.64z~ 2 -0.576z~ 3 

fc 3 =a 3 (3) = -0.576 
BJz) = -0.576 + 0.64 z "' - 0.9z~ 2 + z~ 3 



Mz) = 



A 3 (z)-k 3 B 3 (z) 
\-kl 



= l-0.795z _1 +0.1819z" 



k 2 =a 2 (z) = 0.1819 
fi 2 (z) = 0.1819 



* 2 (z) = 



\-kl 



= l-0.6726z 



^=^(1) = -0.6726 



Ladder coefficient: - 

C 3 (z) = l + 3z~' + 3z~ 2 + z~ 3 
v 3 =C 3 (3) = l 

C 2 (z) = C 3 (z)-v 3 B 3 (z) 

= 1.576 + 236z' 1 +3.9z 2 

C l (z) = C 2 (z)-v 2 B 2 (z) 

= 0.866 + 5.46z~' 
v, = 5.46 

C Q (z) = C l (z)-v l B 1 (z) 

= 4.538 
v = C (0) = 4.538 




3=-0.576 



3= 0.576 



v 3=l 




2=-0. 1819 



2=0.1819 




Fig: Lattice ladder structure. 
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Analysis of sensitivity to quantization of filter coefficient:- 



Let us consider a general FIR filter with system , 



2>. 






H{z) = -^ (1) 



N 



With quantized coefficient coefficient the system function. 

m 

Zb k z~ k 

H(z) = ^^ 



N 



l+2> t z-* 



k=\ 

Where quantized coefficient {b k } can be released to the unquatized coefficients {a k } & {b k } by the 
relation, 

a k =a k + Aa k , k=l,2,3 N (3) 

b k =b k +Ab k , k = 1, 2, 3 m (4) 

Where {Aa k } represents the quantization errors. 

The denominator of H(z) may be expressed to the form , 

D{z) = 1 + ^,1-" =IlU\- p kZ - v ) (4) 

k=\ 

Where, {p k } are the poles of H(z) 

Similarly we can express the denominator of H(z) as, 

D(z) = Ul l (l-p k z- 1 ) (5) 

Where, p k =p k +Ap k , k = \,2 N 

And Ap k is the error or perturtation resulting form the quantization of filter coefficients. 

We shall relate the perturtation errors in the {a k } . 
The perturtation error Ap i can be expressed as, 

Ac-f^ <") 

k=\ (^k 



dp i 8A{z)l Sa k 



5a k {8D{z)l8z) 7= 



(7) 



z=p, 



The total perturtatio error, 



N-k 



a ^=-Zh^ ^* (8) 

The error can be minimized by maximizing the lengths \p t - p\- This can be accomplished by realizing 
the high order filter with either single pole or double pole filter sections. 



Limit cycle oscillation in recursive system: - 

In the realization of digital filter either in digital hardware or software on a digital computer, the 
quantization inherent in the finite precision arithmetic operation render the system non linear. In 
recursive system the non linearities due to finite precision arithmetic operation often cause periodic 
oscillation to occur in o/p even when the input sequence is zero or some non zero constant value. Such 
oscillation in recursive system are called limit cycle and are directly attributable to round off errors in 
multiplication and overflow error in addition. 

Scaling to prevent overflow :- 

In order to limit the amount of non linear distortion it is important to scale the input signal and unit 
sample response between the input and any internal summing mode in the system such that overflow 
becomes a rare event. 

For fixed point arithmetic let us first consider the extra condition that overflow is not permitted at any 
node of the system. Let y k (n) denote the response of the system at k th node when input sequence is x(n) 
and unit sample response between the node and the input in h k (z) . 



y*(«) = 



'^ j h k (m)x(n - m) 



< ^\h k (m)\\x(n-m)\ (1) 



Suppose that x(n) is upper bounded by A k . 
Then, 

|y,(n)|<A, f>,(m)| for all n (2) 

m— — c° 

Now, if the dynamic range of the computer is limited to (-1, 1) the condition, 

\y k (n)\ <l ( 3 ) 

Can be satisfied by requiring that the input x(n) can be such that, 

A,<^ (4) 

For all possible nodes in the system. The condition in (4) is both necessary and sufficient to prevent 
coefficient. 

For FIR filter, (4) become, 

A v <— -^ (5) 

m=0 

Another approach to scaling is to scale the input so that, 

2 

f>,(n)| <C 2 f>(n)| 2 = C 2 £, (6) 

n— —°° n—— °° 

From parsavals theorem, 

£|y»| = — \\H(w)X(wf 



<E t — \\H(wf dw (7) 



From (6) and (7) 



1 " 
C 2 E<E X — \\H(wfdw 



In 
C 2 <— \\H(wfdw (8) 

9 IT J 



-X 

n 



In 



Read the following topics your self:- 

• Representation of Number. 

■ Fixed point representation. 

■ Floating point representation. 

• Rounding and truncation error. 



Design of IIR filter form analog filters: 
Impulse in variance method: 

M A 

tf a (*)=Z— 

<=1 s - Pi 

Taking inverse laplace transform 

M 



Put t = nT 

M 

h a (nT) = Y J Ae P! " T u a (nT) 



i=i 



tf(z) = £/l(")z~" 



= 1 



B=0 

M 



n=0 
M 



-1 

z 



Y^A^uMT) 
i=i 

= 2±A(e p > T z- l ) n u a (nT) 

(=1 n=0 

Af 1 

«=i l-e ft z 



Hence, 



1 1 

>- 



s-p, l-e" ! 'z 



PiT „"1 



Chapter: 6 



Q. Convert the analog filter with the system function H (5) = '- into a digital IIR by means 

(5 + O.I) +9 

of impulse invariance method. 
Solution: 

. + «u . + «u 



(5 + 0.l) 2 +9 (5 + 0.1+ ;3)(5 + 0.1 -y'3) 



1 ■ + - ' 



(5+0.1+7-3) (5 + 0.1-7-3) 

2 2 

1 1 

H (s) = - + - 

(5 + 0.1+J3) (5 + 0.1-J3) 

Then using impulse invariance method, 

1 1 

H(z) = - I- - 



l + e (-0-W3)r z -i 1 _ e (-o. W3 ) rz -, 

l-( e - oir cos3r)z- 1 
\ ] -o.ir OT j -1 . -o.2r T 

1-^ze cos3rjz +e z 

IIR filter design by bilinear transformation: 

IIR filter deign using impulse invariance and approximation of derivative method have a 
sever limitation in that they are approximate only for low pass filter and limited calls of bandpass filter. 

In this we describe a mapping from s-plane to z-plane called bilinear transformation that overcomes 
the limitation of other two design method describe previously. 

Let us consider the analog filter with system function. 

H(s) = — (1) 

5 + a 

Y(s) b 



X(s) s + a 
Or, sY(s) + aY(s) = bX(s) 
Taking inverse laplace transform, we get, 

— — + ay(t) = bx(t) 

dt 

Taking integration on both sides. 
• du(t) 



nT nT 



J dt + a \y{t)dt= \x(t)dt 

m-T nT-T 

Using Trapezoidal rule, 

nT 

\ g(t)dt = T /2[g(nT) + g(nT -T)] 



nT-T 

Then, 



y(nT) - y(nT -T) + aTI l[y(nT) + y(nT - T)] = bT I l[x(nT) + x{nT - T)j\ 
Taking z-transform on both sides, 

Y(z)-z- r Y(z) + aT/2[Y(z) + z~ l Y{z)}= bT / l[x (z) + z~ l X (z)] 
Y(z) __ bT / 2(l + Z- 1 ) 



X(z) 
H(z) = 



1-z' 1 +aTI 2(l+ f 1 



) 



-i\ 



l+z 



+ a 



Hence, s = 



-i\ 



l + z' 



It gives the mapping form s-plane to z-plane. This is called bilinear transformation. 
To investigate the characteristics of bilinear transformation let, 



, ./"■ 



Now 



z = re J 

S = (T+ jQ 

2fz-0 



s = ■ 



z + l 



re 



./"■ 



re 



+ 1 



( '2 



r -l+j'2rsinw 
r + 2r cos w + 1 



Comparing with, 
s = cr+ jQ. 

2 ' 
o = — 
T 

2 



r 2 +2rcosw + l 



a = 



2r sin w 

v r 2 + 2rcosw + l 



If r <1, then a < and if r >1, then, a > consequently, the LHP is S maps into inside a unit circle in 
z-plane and RHP in s maps into outside a unit circle. 
When, r= 1, cr = 0, 

2( 2sinw ' 



Q = 



v 



2 + 2cosw 



Q. = — tan(w/2) 



w = 2 tan 



Q.T 



Itan-'toT^) 




Fig: mapping between frequency variables w and Q. resulting from bilinear transformation. 

# Convert the analog filter with system H (s) = ^ into digital IIR filter by means of 

(5 + O.I) + 16 

bilinear transformation the digital is to have resonant frequency of w r = n 12. 

Solution: 

The resonant frequency of analog filter is Cl r = 4 . & resonant frequency of digital filter is 



Q. r = — tanw r 12 

. 2 n 

4 = — tan — 

T 4 



T = 



1 



fl-z^ 
[l+z- 1 , 


= 4 


yl+Z- 1 } 



Now, s = 2/T 



Using bilinear, the system function of digital filter becomes H(z) = 



0.125 + 0.006z _1 -0.118z" 2 



\-z 

l + z~ 



+ 0.1 



1-z 

1+z" 



+ 0.1 



+ 16 



H(z) = 



l + 0.0006z _1 +0.95z _1 
0.125 + 0.006z _1 -0.118z" 2 



l + 0.95z 



Q. Design a single pole low pass digital filter factor with 3 dB bandwidth of 0.2 n using the bilinear 

Q 

transformation apply to the analog filter H(s) = — . Where, Q. c is a 3dB bandwidth analog filter. 



s + Q, 



Solution: 

W c = 0.27T , 

For analog filter. 

2 
Q, c = — tanw c /2 

= |tan(0.1Q( )) 
0.65 



Analog filter has system function 
H(s ) = °- 65/T ' 




H(z) = 



2.65-1.35z _1 

0.245(1 + z') 
l-0.509z ' 



Matched z-transformation:- 

Another method for converting analog filter into equivalent digital filter is to map the poles and zeroes 
of H(s) directly into poles and zero in z-plane. 

Suppose the transfer function of analog filter is expressed in the factored formed. 

U k=l {s-p k ) 

Where, \z k } and {p k } are zeroes and poles of the analog filter. 
Then system function for digital filter is 

rCa-^V 1 ) 



H(s) = 



nLa-^V 1 ) 



Thus each factor of the form (s-a) in H(s) is mapped into the factor (1 - e" z ) . This mapping is called 
matched z-transformation. 



# Convert the analog filter with system function H a (s) = ' — - into digital IIR filter by matched 



(s + O.iy+9 



z-transformation method. 

5 + 0.1 



(5 + 0.1) 2 -(3j) 2 ) (5 + 0.1 + 37)(5 + 0.1 -;3) 
Using matched z-transformation method. 

') 



H(z) 



¥■ 



L -o.it -1 
\\--e z 



. e (-0.i-3 J )r z -iV 1 _ g (-0.i +J -3)r z -, 



1 



/! o.ir -K 



-e e J z )(l-e e z ) 



O.IT i"3T -1- 
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Butter worth filter: 

I |2 1 

We have, 7"„(jw) = 



l + £-V" 



When £ 2 =1 

It t ■ \\ 2 * 

1 + 6) 

This function is known as Butterworth response. From this equation we observe some interesting 
properties of Butterworth response, 

(1) The Butterworth filter is an all -pole filter. It has zero at infinity (a> — > °°). 

(2) |r„(jO)| = lforalln. 
1 



(3) \T n (jl)\ = —j= = 0.707 for all n. corresponding to 3 dB. 

V2 

(4) For large co, \T n {jaJ)\ exhibits n-pole roll-off. 

|Tn(jca 




Fig: Maximally flat response. 



Butterwoth filter : - 

a(G)) = -20\og\T (jw)\dB 

Low pass filter specifications, 
Passband frequency = Q. p 

Stopband frequency = Q s 
Passband attenuation = a p 
Stopband attenuation = cx s 
Order of filter, 

n = /* — ; — r 

2iog(n,/nJ 

3-dB cutoff frequency, 
Q 



2Q 01 »p _1 )2n 



Chebyshev filter: 

We have chebyshev magnitude response, 



Kun = 



i 



l + £ 2 C 2 n (co) 



(1) 



C n {co) = cos(«cos l co) for \co<l\ 

= cosh(ncosh _1 fy) for |<y| > 1 (2) 

With equation (2), the transfer function magnitude in (1) is determined for all values of co . The function 
is plotted for n = 6. 



n odd begin here 



n even begin here 




Fig: Sixth order chebyshev Tp magnitude 



(1) C„(0) = for odd and C„(0) = 1 for n even. 
|r n (jO)| = l for n odd. 
1 



KU0) = 



vr 



for n even. 



+ £ z 



(2) At m = \, C„(l) = l 



il 



for all n. 



+ £ 

f 1 nO.la, 1 A 



cosh 



10" 



1 



n = 



v 10° lff "-l y 



cosh"'(Q s /^J 

equal ripple filter. 

High attenuation in stop band and steeper roll-off near the cut-off frequency. 



* For Q, p = 1, Q s = 3.33 and a p = 03dB a s = 22dB Which filter is best. 

Solution, 

For butterworth filter n = 5. for chebyshev filter n = 3. 
Greater efficiency of chebyshev filter compared with Butterworth filter. 



Frequency transformation :- 

If we wish to design a high pass or bandpass or bandstop filter it is a simple method to take a low pass 
prototype filter (butterwoth , chebyshev) perform a frequency transformation. 

One possibility is to perform the frequency transformation in analog domain and then to convert 
analog filter into corresponding digital filter by a mapping of s-plane into z-plane and alternative 
approach is first to convert the analog low pass filter into digital low pass filter into a desired digital 
filter by a digital transformation. 

In general, these two approaches yields different results except for bilinear transformation in which 
case the resulting filter designs are identical. 



Frequency transformation in analog domain: 



Type of transformation 



Low pass 



High pass 



Band pass 



Band stop 



Transformation 



5^ 



^ 

a 



S ^5A 



s^a. 



s 2 + Cl L Cl u 



s^a. 



g(fl M -fl L ) 



Band edge freq. of new filter. 



«,■ 



«,■ 



«L A 



«L A 



Q. Transform the single pole low pass butterworth filter with the system function H(s) = 
band pass filter with upper and lower band edge frequency Q. u and Q L respectively. 



£1. 



s + Q.. 



into 



Solution: 



s = n. 



H(s) = 



* 2 + ^A 






*("u-"J 



5 2 + (Q„-Qj5 + Q L Q„ 

Frequency transformation in digital domain (spectra transform): 

Conversion of prototype digital low pass filter into either band pass, band stop, high pass, 
or another low pass digital filter. 

Need of spectral transformation :- 

Impulse invariance method and mapping of derivative are inappropriate to use in designing 
high pass and many band pass filter due to a aliasing problem. 

Consequently one could not employ an analog frequency transformation followed by 
conversion of the result into digital domain by use of this two mapping instead it is much better to 
perform the mapping analog low pass filter into digital low pass filter by either of these mapping and 
then to perform frequency transformation in digital domain. 
The problem of aliasing is avoided. 



Type of transformation 



Transformation 



parameter 



Low pass 



_i z ' - a 
z -> T 



CO = band edge frequency of new filter 



l-az 



sin 



a = 



a?-® pi 



sin 



r G> p +G>S 



High pass 



z -> — 



z l + a 
l + az' 1 



cos 



a = -- 



r co p +o) p ^ 



COS 



f • \ 

co p -co p 



Band pass 



-2 -1 , 

_i z —a,z +a 7 

Z -> y^ p^ 

a 2 z -a x z +1 



co L = Lower band edge frequency 

CO u = Upper band edge frequency 

2ak 

a, = 

k + l 

k-1 

a = 

k + l 



a = 



k = cot 



cos 




r co u +co L ^ 




V 2 j 


COS 




®u-®L 




V 2 , 






.tan 



k/2) 



y 



Band stop 



_i z 2 -a,z 1 + a, 
z -> 



"i" 



a 2 z -a x z +1 



2a _ 1 - & 

, a 2 — 



fc + 1 



1 + Jfc 



a = 



k = tan 



cos 




r o) u +co L ^ 




V 2 , 


COS 




f ' A 




V 2 , 



CO.. - CO, 



tan 



k/2) 



# Design a low pass butterworth filter to meet the following specifications: 
Passband gain = 0.89 
Passband frequency = 30 hz. 
Attenuation = 0.20 
Stopband edge = 75 Hz. 

Solution: 
Q p = 2^x30 = 188.4 rad/sec 

Q. s =2^x75 = 471.2 rad/sec 

a p =-20\og(0.89) = l.0ldB 

a s = -201og(0.20) = U,9SdB 



log 



n = 



^1Q°- Ig --1 N 
v 10 aia '-l, 



= 2.466 « 3 



21og(Q p /Qj 
3 -dB cut off frequency, 

". =7 — ^ \ = 23.55 rad/sec 

(lO 01a "-l) 

From table for N = 3, and Q c =1, we have 

#00 = 7 \ t 

(s + l)/(s 2 + s + l) 



This function is normalized for Q c = 1. However Q c = 235.55, we need to denormalized H(s) by 
Q. =235.55 rad/sec. 



H(s) 



1 



s = s/ci L { s /a c +i)(s 2 /n 2 c + s/a c + i) 

Where, Q c = 235.55 rad/sec. 



# Design a low pass FIR filter with specification co p = 0.2k , 0) s = 0.65;r, a p = OAdB , a s =l.5dB . Use 

bilinear transformation method. 
Solution: 

a p = -t<m(a) p /2) = -tan(0.2;r/2) = 0.649(7 = 1)) 

Q,=-tan(ft> s /2) = 3.263 
Order of filter, 
log 

2\og(Cl p /Cl s ) 
3 -dB cut off frequency, 
CI. 



n = 



^ io oig --r 

10 o.i«„ _ l 
v U ^-1.783-2 



& =■ 



10°- Ia '-lp» 



= 1.164 



From table for n = 2 and £l c = 1 , we have 

H(s) = - 2 

s 2 + 1.4145 + l 

Now denormalizing with Q c =1.164 

II(s) ] 



s = s/Q s 2 /Cl 2 +lAls/Cl+l 



a 2 



s 2 + lAUsQ c + Q 2 
Q 



s 2 + 1.414 sQ_ + Q 



Now using bilinear transformation, 



s = ■ 



T 
H(z) = 



U + ly 


= 2 


vz + U 



1.354 



fz-l\ 



+ 1.645 



,^1 



+ 1.3548 



Q. Convert a given H (5) = 



0.2 



5 + 0.2 

with pass band edge frequency Cl p = 0.5 
Solution: 

0„Q„ 0.2x0.5 _ 0.1 

s s 

0.2 0.25 



of band edge frequency Cl p = 0.2 into a high pass filter H n (s) 



S^ 



H n (s) = H, 



'n,<C 



0.1/5 + 0.2 0.25 + 0.1 



*W = 



5 + 0.5 



Q. Use bilinear transformation to obtained digital low pass filter to approximate H(s) = 

Assume cut off frequency of 100 Hz and sampling frequency of 1 khz. 
Solution: 

2^-xlOO 



5 2 + V25 + 1 



G)„ = 



(Normalizing with f s ) 



1000 

= 0.2;r rad/sec 

n c =|tan(fl> c /2) = 0.65(r = l) 

Now, denormalizing H(s) with Q. c = 0.65 we get, 

1 

5 = 5/Q 



H L (s) = H(s) 



\ 



0.65 



+ V2(5/0.65) + l 
0.4225 



5 2 + 0.9195 + 0.4225 



Now we substitute, 5 = 



H(s) 



z-l 

z + K 

0.4225 



= 2! 



z-l 

z + l 



(r = i) 



2 


vz + U 


-2 

+ 0.919 


2 


'z-lY 



+ 0.4225 



Q. Using bilinear transformation design a butterwoth filter which satisfy the following condition. 



0.8 < H(e jw ) <1 

\H(e iw )\<0.2 
Solution: 



0<a;<0.25 
0.6 < a <n 



|Tn(jco)| 



1 

0.8 
0.2 



a p = -201og(0.8) = 0.93dB 
a s =-20\og(0.2) = 13. 97 dB 
co p = 0.2k , co s = 0.6;r 

a„ =— tmico / 2) =0.65 

P rji \ p / 

2 / \ 
a s =—tm{co s /2) = 0J5 



log 



10° 



1 



,0.1ar„ 



n = 



Q. = 



10"^ -1 

2io g (n p /nJ 
ft. 



= 0.75 



(io olff "-i) 2 " 

From table, for n = 2 and Q. c = 1 

H(s) = - 2 1 

s 2 + 1.414s + l 

Now, denormalizing with Cl c = 0.75 

H(s) l 



Transition band 




s = s/Q c (s/0.75) 2 + 1.414(s/0.75) + l 
0.56 



s 2 + 1.065 + 0.56 
Using bilinear transformation, 

2fz-l A 



5 = 



H(s) = 



z + 1 



z-1 
.z + 1 



(r = l)) 

0.56 

-2 

+ 1.06 2 



z-1 

z + 1 



+ 0.56 



Stop band 



C0 S = 0.27t ms = 0.6jc 



CD 



Q. Convert low pass butterworth filter with system function H(s) = — r- into bandpass filter 

l-0.509z _1 

with upper and lower cut off frequency 03 u and 03 c respectively. The low pass filter has 3 dB bandwidth 

Q) p = Q.Ijc. [o3 u = 3,7V I '5, 0) L = 2/r/5) 

Solution: 



i G>„ ~ 03, 

k = cot " L 



= 1 



tan(Q) p /2) 



cc = 



COS 


~03 u 


+ 03 L ~ 




2 


COS 


03 u 


-03 L 




2 



= 



-2« t 
k + \ 
k-\ 
k + l 



= 



= 



Now, z'^-z 2 



TF becomes, 



-2> 



^(z) = 



0.245(1 -z' z ) 
l + 0.509z" 2 
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Chapter: 7 

FIR filter design:- 

FIR filter:- The digital filters can be classified either as finite duration unit impulse response (FIR) 
filters or infinite duration unit impulse response (IIR) filters depending upon the form of unit impulse 
response of the system. In FIR system the impulse response sequence is of finite duration. 
For example the system with impulse response, 

1 2 for \n\ < is FIR system 
h{z) = \ 

[0 otherwise 

FIR filters are generally implemented using structures with no feedback (i.e non recursive structure). 
An FIR filters of length 'N' can be describe by the following difference equation. 
y(ri) = b x(n) + b x x(n - 1) + b n _ x x{n -m + l) 

M-\ 

or, y(n) = ^b k x(n - k) Where, bk is filter coefficient. 

k=Q 



Symmetric and antisymmetric FIR filters: 

The unit sample response of FIR filter is symmetric if it satisfies the following condition. 
h(n) = h(M -l-n) , n = 0, 1,2, M-l. 



As in example let us consider the unit sample response given in figure for M = 6 sample. 

h(n) 



6 6 

9 9 



12 3 4 5 



h(n) = h (5-n) 
h(0) = h(5) =2 
h(l) = h(4) = 4 
h(2) = h(3) = 6 



This is symmetric FIR filter. 

The unit sample response of FIR filter is antisymmetic if it satisfies the following condition. 

h(n) = -h (m-l-n) , n = 0, 1, M-l 

As an example let us consider the unit sample response shown in figure for M = 6 samples. 

h(n) 



1 2 



4 5 



-4 



h(n) = -h(5-n) 

h(0)=-h(5) =2 

h(l)=-h(4)=4 

h(2)=-h(3)=6 

This is antisymmetric FIR filter. 



Note: The condition for linearity phase FIR filters is h(n) = ±h(M -\-n) 



The phase of linear phase filter is e 



-J<»\ 



M-\ 



# Show that the digital FIR filter with impulse response h(n). 

h{n) = {2,4,6,6,4,2} is liner phase system. Is this antisymmetric? 



Solution: 

h(n) = ±h(M -l-n) 
Here, M = 6. 
h(n) = ±h(5-n) 

Where, 
h(0) = 2, h(l) = 4, h(2) = 6, h(3) = 6 
h (4) = 4, h(5) = 2, 
h(0) = h(5) = 2 
h(l) = h(4) = 4 
h(2) = h(3) = 6 

Hence the system has linear phase and h(n) is not antisymmetric since h(n) = h(M-l-n) 

# A digital filter has impulse response given by h(n) = {1,0, 0, 0, 0, 0, -1}. What is its system function. 
Which class of linear phase filter does this system belong to ? Justify. 

Solution: 

System function, 

H(z) = f j h(n) Z - n 

= l-z 

Linear phase FIR filter satisfies the conditions, 

h(n) = ±h(M -l-n) 
Here, M= 7, h(n) = ± h(6-n) 

h(0) = -h(6) = 1 

h(l) = -h(5) = 

h(2) = -h(4) = 

h(3) = -h(3) = 
Hence, this system belong to antisymmetric. 

# A linear phase filter has a phase function e~' w . What is the order of the filter. 

Solution: 

The phase of linear phase filter is given by, 



Comparing with e 
M -1 



e 

j2w 



' M-\ 



= 2-->M =5 

2 

Design of FIR filter using fourier series method: - 

The frequency response H d (e JW ) of a system is periodic in In .We know that , periodic function can be 
expressed as linear combination of complex exponentials. Therefore the desire frequency response of FIR filter 
may be represented by fourier series. 

H d (e ja) )= f>(n)<T M (1) 

n— —oo 

Where, the fourier coefficients h d (n) are the impulse response sequence of the filter given by, 

h d {n) = \* H d {e im )e im dco (2) 

Also, z-transform of the sequence is given by, 



H(z)= J>»Z- 



(3) 



Equation (3) represents a non-causal digital filter of infinite duration to get an FIR filter transfer function the 
series can be truncated by assigning. 

h A yn) for \n\< 

h{n) = \ dK J J ll 2 (4) 

otherwise 

Then we have, 

H(z)= j^H^z- 1 



.(5) 



h 



^ 2 j 

N-\ 

2 



JV— 1 



+ h(l)z~ l + h(0) + + h 



r N -\\ 



\ ^ j 



(JV-V2) 



H(z) = h(0)+Y j [h(n) Z - n + h(-n)z"] 



(6) 



For symmetrical impulse response symmetrical at n = 0. 

h(-n) = h(n) (7) 

Therefore equation (6) becomes, 

N-\ 

H(z) = h(0)+^h(n)[z" + z- n ] (8) 

n=\ 



This TF is not physically realizable. Realizability can be brought by multiplying the equation (8) by z 

N-\ 
Where, is delay in samples. 



N-\ 



We have, 



H\z) = z 



N--H (z) 



= z 



- lN -2 



JV-1 

2 



h(0) + ^h(n)[z" + z~"] 



(9) 



# Design ideal lowpass filter with following frequency response. 

\\ for-nl2<co<nl2 

|p forKll<(0<K 
Determine the values of h(n) for N = 11. Find frequency response. 



H d (e ] ) 



Solution: 



W = T-f H d (e m )e im dco 

2k j -* 



1 f!tl2 

= — \ !■ 



jam 



dco 



HdCeiw) 



-tt/2 



ji/2 



sinn — 

2_ 

tin 



- °° < n < oo 



Now truncating h d (w) to 1 1 samples, we get, 

-7T 



n(n) = 



sinn 



— forln < 5 

;zn 

otherwose 



h(0) = '/2 

h(l) = h(-l) = 0.3183 
h(2) = h(-2) =0 
h(3) = h(-3) =-0.106 
h(4) = h(-4) =0 
h(5) = h(-5) = 0.06366 



N-\ 

2 



TF, H(z) = h(0) + 2>(n)[z" + Z~"] 

«=i 

= 0.5 + 2n(n)[z" + z""] 

= 0.5 + 0.3 1 83CZ 1 + z" 1 ) - 0. 106(z 3 + z~ 3 ) + 0.06366(z 5 + z~ 5 ) 
TF, of realizable filter will be, 

H\z) = z [ 2) H(z) = z- 5 H(z) 

= 0.06366 - 0.106z 2 + 0.318z~ 4 + 0.5z~ 5 + 0.318z~ 6 - 0.106z~ 8 + 0.06366z 10 
From above, the filter coefficients of causal filer are given as, 
h(0) = h(10) = 0.06366 
h(l) = h(9) = 
h(2) = h(8) = -0.106 
h(3) = h(7) = 
h(4) = h(6) = 0.013183 
h(5) = 0.5 

Window techniques for designing FIR filters: 

In this method we begin with the design frequency response specification H d (w) and determine the 
corresponding unit sample H d (n) indeed H d (n) is related to H d (w) by the fourier transform relation. 



H d (QJ) = Y,h d (n)e~ im (D 

n=() 

Where, 

h d {n) = ^-\ K H d {co)e im dco (2) 

In general unit sample response H d (n) obtained form equation (2) is infinite in duration and must be truncked at 
some point say at n = M-l, to yield FIR filter of length 'L'. Truncation of H d (n) to a length M-l is equivalent to 
multiplying H d (n) rectangular window defined as, 

[\ n = 0,1,2 M-l 

co(w) = \ (3) 

[0 otherwise 

This unit sample response of FIR filter becomes , 

[h, n = 0,l M-l 

h(n) =h d (n)co(n) = \" , . (4) 

[0 otherwise 

It is instructive to consider the effect of window function on the desire frequency response H d (CO) . The 

multiplication of the window function C0(w) with h d (n) is equivalent to convolution of H d (CO) with W(co) . 

Where W(0)) is the frequency domain representation of window function. 

M-l 

W(60) = ^W(n)e- jm (5) 

Thus the combination of H d ( CO ) with W(co) yields the frequency response of FIR filter. 

H(co) = — I" H d (v)W(co-v)dv 

IK*-* 



rier transform of rectangular window is , 


W(co) 


M-l 
X ' - icon 

_\-e- jaM 
\-e-' m 





-Mm -i/2) sin(coM/2) 

sin(co/2) 

The window function has a magnitude response, 

, | |sin(ftM/2| 

WC^)^ 1 ^ -, — -hrr 1 -K<C0<7Z 

fva\coj2\ 
And a piece wise linaer phase . 

e(co) = oi — —1 When, sin(ffiM/2)>0 

I 2 J 

= -c4 +K When, sin(ft»M/2)<0 
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Name of window 



Time domain sequence 

ffl(n), 0<n<-l 



Rectangular 



1 



Hamming 



0.54-0.46cos- 



2m 
M-l 



Hanning 



If, 2m 

1-cos- 



M-l 



Kaiser 



A M-lY ( M-0 



v 



a 



M-l 



Where I is zero order Bessel function. 



J 

^ 



•rf" 



Rectangular 



• \ Hamming 
\ Hanning 



M-l 



kaiser 



M-l 



Name of window 



Window function. 



0){n), \n\ < 



v 



j 



Hamming 



0.54-0.46cos 



27m 
M-l 



Hanning 



If, 2m ^ 

—\ 1 + cos 

2l M-l 



Gibbs phenomenon in FIR filter design: 




Fig: frequency response of filter. 



frequency 



We observe that relatively large oscillation or ripple occur near the band edge of filter the oscillation increases in 
the frequency as N increases but they do not diminish in amplitude. Large oscillation are direct result of large side 
lobes existing in frequency characteristics of W(co) of rectangular window. Fourier series representation of 

H d (CO) , multiplication of h d (CO) with rectangular window is identical to trunked in fourier series representation 
of desire filter characteristic H d {O0) . The truncation of fourier series is known to introduce in frequency response 
characteristic H(oo) due to non uniform convergence of fourier series at discontinuity. The oscillatory behavior 
near the band edge of the filer is called Gibbs phenomenon. 

To alleviate the presence of large oscillation in both the pass band and stop band we should use a window 
function that contains a taper and decay towards zero gradually instead of abruptly. 

Design of FIR filter by Frequency sampling method: - 

In the frequency sampling method of IIR filter design we specified the desire frequency response H d (CO) at a 
set of equally spaced frequency. 

2n 

Namely CO, = (k + a) (1) 

M 

Where, k = 0, 1, — for M odd. 

2 

M 

K = 0, 1 1 for M even. 

2 

a = 0, or — 

2 

And solve the unit sample respose n/w of the FIR filter from these equally spaced frequency specification. Now 
the desire frequency response is 

M-l 

H(w)= ^h(n)e~ im 

Suppose that we specify the frequency response of the filter at the frequency given by equation (i). 

'2n N 

[k + a) 

j 



H(k + a) = H\ — (k + a) 
M 



H{k + a)=Y J h(n)e- i2Mk+a)n ' M (2) 

Where,k = 0, 1, M-l 

It is simple matter to invert (2) and express h(C0) in terms of H(k + CX) 

i M-l 

h(n) = — Y J H(k + a)e j2 ' r(k+a)n ' M (3) 

^" k-o 

Where, n = 0,1 M-l 

This relationship in (3) allows us to compute the values of unit sample response h(n) from the specification of 

frequency sample H(k + a) , k = 0, 1, M-l. Note that when CX = (2) reduces to discrete fourier transform 

(DFT) of the sequence h(n) and the expression (3) reduces to IDFT. 

# Design a low pass FIR filter with 1 1 coefficient for the following specification. 

Passband frequency = 0.25 khz 

Sampling frequency = 1 khz. 
Use rectangular window, Hamming window and hanning window. 
Solution: 



f c = 0.25khz 



co = 



2;rx0.5 n 



1 (Till 

h,{n) = — f l.e jm dQ) 



1 sin(nff/2) 

2 {jut 1 7) 



1 



H(ro) 



-7t/2 n /2 tj 



hd(0) = Vi 

h d (l) = h d (-l) = 0.3148 
h d (2) = h d (-2) = 
hd(3) = hd(-3) = -0.0162 
h d (4) = h d (-4) = 
h d (5) = h d (-5) = 0.06369 

Rectangular window :- 

w{n) = 1 - 5 < n < 5 

h(n) = h d (n) w(n) 
h(0) = h d (0) w(n) = 0.5 
h(l) = hd(-l) =0.3184 
h(2) = h d (-2) =0 
h(3) = h d (-3) =-0.1062 
h(4) = h d (-4) =0 
h(5) = hd(-5) =0.06369 



Hamming window: 

W ham (n) = 0.54+0.46 cos 



/ 



= 0.54+0.46cos 



2701 

M-\ 

r 2m^ 



\ 



\- 



j 



v 



10 



j 



h{n) = h d {n)W ham {n) 
W ham {0) = l 

W ham {D=W ham {-l) = 0.9121 
WL.(2) = W r ta ,(-2) = 0.6828 
WL,(3) = W tall (-3) = 0.3970 
WL.(4) = W ham {-4) = 0.1679 
^ flm (5) = W /Mffl (-5) = 0.08 



<n< 



J 



V 2 
5<;r<5 



h(0) = h d (0)W ham (0) = 0.5 

h{\) = h{-\) = 0.2904 

h{2) = h(-2) = 

h(3) = h(-3) = -0.04225 

h(4) = h(-4) = 

h(5) = h(-5) = 0.005095 

H(z) 

Similarly Hanning window.. 



Q. A low pass filter is required to be design with design frequency response which is expressed as follows, as we 

\e' i2(0 -nl2<\co\<nl\ 

[0 for x/4<\a\< n 

Obtain the filter coefficient h d (n) if the window function is defined as, 

[\ < n < 4 
W(n) = { 

[0 otherwise 

Q. Design a low pass filter having desire frequency response given as, 

\e-»°> Q<co<7tl2 

[0 7tl<ca<7t 

Obtain filter coefficient h(n) for M = 7 using frequency sampling method. 
Solution: 

W.=—(k + a) 

k M 

a = 

1 6 
h(n) = — y_ i e J e J n = 0, 1, 6. 

' k=Q 
' k=Q 



# Remez exchange algorithm (prokish book «IIC 3TTqT H(i<H) 



Reference for chapter (8) 

TMS 320 Texas instrument (Sanjaya sharma) 

Bit serial arithmetic (Digital logic book) 

Serial adder. 

Distributed arithmetic and pipelined implementation. 



